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Abstract 


A  co-axial  contra-rotating  helicopter  in  forward  flight  has  stall  mechanisms  that 
are  intrinsically  different  from  those  on  a  traditional  helicopter.  Traditional  helicopters 
need  a  cyclic  pitch  mechanism  to  balance  the  rotor  lift  about  each  rotation,  which  leads  to 
dynamic  stall  from  the  rapid  oscillations  in  pitch.  Co-axial  contra-rotating  helicopters, 
which  have  a  fixed  pitch  about  each  cycle,  encounter  a  sinusoidal  oscillation  in  Mach 
number  with  the  mean  velocity  seen  as  the  rotational  velocity  while  the  half  amplitude  is 
equal  to  the  forward  flight  speed.  The  fluid  dynamic  mechanism  limiting  the  forward 
flight  speed  is  entirely  different  from  that  of  traditional  dynamic  stall  studies.  This  work 
sets  out  to  design,  create,  and  study  the  application  of  an  oscillatory  compressible  flow 
field  on  a  NACA  0012  airfoil  in  order  to  experimentally  model  this  flow  situation.  A  few 
different  pressure-sensitive  paints  and  imaging  techniques  were  developed  for 
investigation  of  this  oscillatory  effect.  Ultimately  a  fast-acting  bi-luminophore  pressure- 
and  temperature-sensitive  paint  was  chosen  for  the  investigation  which  uses  a  polymer- 
ceramic  basecoat  and  a  mixture  of  luminescent  elements.  These  measurements  were 
made  using  a  two-camera,  single-shot,  intensity-based  pressure-sensitive  paint  technique. 
Temperature-corrected  pressure  measurements  were  made  and  accounted  for  the  intrinsic 
temperature  sensitivity  of  pressure-sensitive  paint.  This  dual-luminophore  technique 
allows  for  accurate  unsteady  pressure  measurements  in  a  non-uniform  and  varying 


temperature  environment;  however,  due  to  inadequate  unsteady  surface  pressure  tap 
measurements  the  pressure  results  were  limited  to  the  steady  runs. 

This  work  involved  the  design  and  creation  of  a  modification  to  Ohio  State’s  6”  x 
22”  Transonic  Wind  Tunnel  to  enable  oscillations  of  the  freestream  Mach  number.  The 
test  section  Mach  number  in  this  blowdown  facility  is  dynamically  set  by  varying  the 
choke  area  downstream  of  the  test  section  in  an  oscillatory  fashion.  The  current 
configuration  of  the  facility  can  produce  Mach  number  oscillations  between  0.44  and 
0.64  for  a  Reynolds  number  range  of  17  -  43  million  per  meter  at  oscillation  frequencies 
up  to  2 1  Hz.  Unsteady  measurements  in  this  oscillating  freestream  flow  have  been  made 
on  a  NACA  0012  airfoil. 

Unsteady  shock  location  measurements  were  made  at  angles  of  attack  of  9,  10  and 
11  degrees  and  frequencies  of  2.1,  9.5,  15.25  and  21  Hz  on  the  NACA  0012  airfoil. 
Detailed  measurements  of  the  shock  movement  were  made  with  these  advanced 
measurement  techniques  in  order  to  investigate  unsteady  effects  of  the  oscillatory 
freestream  flow.  Unsteady  effects  were  pronounced  for  a  reduced  frequency  of  0.037 
which  is  below  the  typical  quasi  steady  to  unsteady  threshold  of  0.05.  Coefficient  of 
pressure  measurements  for  the  steady  runs  were  validated  with  historical  data.  It  was 
found  that  the  coefficient  of  lift  measurements  were  in  very  good  agreement,  while  the 
coefficient  moment  had  significant  errors.  Furthennore,  the  PSP  measurements  were 
compared  with  particle  image  velocimetry  data  of  other  researchers  in  order  to  form  a 
comparison  of  the  on-  and  off-body  fluid  dynamics  at  a  frequency  of  9.5  Hz  and  angles  of 
attack  of  9  and  10  degrees.  The  single-shot  pressure-sensitive  paint  technique  was  able  to 

iii 


measure  buffeting,  which  was  found  to  be  highly  three  dimensional  over  the  span  of  the 
airfoil.  Similar  shock  location  unsteadiness  due  to  buffeting  was  also  measured  in  the 
forced  oscillation  cases  at  lower  Mach  numbers  than  steady  runs,  causing  aperiodic 
behavior  at  certain  azimuth  locations  with  higher  Mach  numbers. 
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Chapter  1:  Introduction 


Section  1.1:  Helicopters  in  Forward  Flight 

When  a  traditional  helicopter  is  in  forward  flight  the  rotor  blade  experiences  a 
difference  in  relative  velocities  between  the  advancing  and  retreating  blades,  as  illustrated 
in  Figure  1.  Therefore,  a  cyclic  rotor  pitch  adjustment  is  imposed  to  balance  the  lift 
distribution  about  the  rotor  azimuth  {¥),  such  that  the  angle  of  attack  is  low  on  the 
advancing  blade  and  high  on  the  retreating  blade.  Traditional  helicopters  are  limited  in 
their  forward  flight  speed  by  dynamic  stall  or  shock-induced  stall.  Dynamic  stall  is 
caused  by  the  rapidly  changing  angle  of  attack  (a)  between  the  advancing  and  retreating 
rotor  blades  from  the  need  to  balance  the  lift  distribution.  During  high-speed  flight, 
maneuver,  and  at  high  altitude,  this  rapid  pitching  motion  to  post  stall  angles  of  attack 
and  oscillatory  relative  flow  field  causes  dynamic  stall.  Figure  2  illustrates  a  traditional 
UH-60  Blackhawk  helicopter  in  forward  flight. 
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Figure  1 :  Schematic  diagram  of  a  helicopter  in  forward  flight,  illustrating  the  relative 
velocity  of  the  rotor  blade  as  a  function  of  its  azimuth  position 


Due  to  this  dynamic  stall  limitation,  non-traditional  helicopter  configurations  have 
arisen  from  the  need  to  fly  faster  or  more  efficiently.  One  method  to  avoid  the  asymmetry 
of  lift  is  the  use  of  contra-rotating  blades;  examples  include  the  Kamov  Ka-52,  and  more 
recently  the  Sikorsky  X2,  as  seen  in  Figure  3  and  Figure  4  respectively.  With  these 
contra-rotating  configurations  it  is  no  longer  necessary  to  have  a  tail  rotor.  Additionally, 
cyclic  pitch  mechanisms  are  not  necessary  to  balance  the  aircraft’s  lift  about  the  rotor 
azimuth  for  co-axial  helicopters.  Therefore,  the  rotor  angle  of  attack  is  fixed  throughout 
its  revolution. 
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The  only  phenomena  limiting  the  forward  flight  speed  of  coaxial  contra-rotating 
helicopters  are  shock-induced  stall,  power  available,  vibratory  loading  and  the 
aerodynamic  drag  performance  of  the  helicopter.  In  order  to  engineer  rotor  designs  that 
expand  the  operating  envelope  of  these  rotorcraft,  it  is  imperative  to  investigate  the  flow 
mechanisms  that  lead  to  shock-induced  stall  on  an  airfoil  subject  to  a  time-varying 
compressibility  condition.  This  leads  to  a  need  for  advanced  measurement  techniques 
which  have  the  ability  to  track  the  motion  of  standing  shock  waves  with  greater 
resolution  than  traditional  pressure  taps.  In  particular,  there  is  a  limited  amount  of 
knowledge  about  the  impact  of  unsteady  oscillatory  compressible  flow  effects  on  shock- 
induced  stall  that  is  typical  of  contra-rotating  helicopter  operations  in  maneuver  or  high¬ 
speed  forward  flight. 
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Figure  2.  UFI-60  Blackhawk  in  forward  flight1 
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Figure  3:  Kamov  Ka-52  coaxial  contra-rotating  helicopter 


Figure  4:  Sikorsky  X2  compound  coaxial  contra-rotating  helicopter 
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Section  1.2:  Fluid  Dynamic  Background 


Numerous  investigations  have  focused  on  the  mechanisms  of  dynamic  stall,  for 

4  8* 

example  McCroskey  et  ah,  McAlister  et  al.  and  Tijdeman  and  Seebass;.  '  Typically  these 
studies  have  been  performed  in  the  low-speed,  incompressible  regime  (M<0.1).9'"  These 
investigations  have  focused  on  the  leading-edge  vortex  that  is  shed  from  the  airfoil  which 
leads  to  lift  and  moment  dynamic  stall,  while  employing  mechanisms  that  vary  either 
airfoil  pitch  or  freestream  velocity.  Saxena  utilized  hot  wires,  surface  pressure 
measurements  and  surface  hot  film  gages  along  with  silk  tufts  to  conduct  a 
comprehensive  on-  and  off-body  study  of  a  NACA  0012  airfoil  in  an  oscillating 
airstream.  ~  These  measurement  techniques  were  used  to  investigate  relatively  high 
reduced  frequencies  (k=)  of  0.18  to  0.9  and  a  chord  Reynolds  number  of  250,000.  The 
thrust  of  the  investigation  was  on  the  effects  of  angle  of  attack  on  the  steadiness  of  the 
flow  field,  concluding  that  the  flow  was  quasi  steady  above  and  below  the  critical  (or 
stall)  angle  of  attack.  However  it  was  found  that  near  the  onset  of  stall  the  airfoil  was 
more  prone  to  perturbances  and  oscillated  between  stalled  and  attached  flow  over  the 
airfoil  with  the  flow  oscillation.  At  angles  below  the  static  stall  angle,  trip  strips  had  no 
effect  on  the  outcome  and  it  was  concluded  that  Reynolds  number  and  transition  effects 
had  no  effect  on  the  flow  phenomenon. 

Selerowicz,  Szumowski  and  Meier  in  two  separate  studies  investigated  forced 
background  flow  oscillations,  0.65  <  M  <  0.75,  at  the  buffeting  frequency  and  angle  of 
attack  of  a  NACA  0012  airfoil  using  interferometry  and  pressure  taps.13’ 14  However,  with 
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only  14  Kulite  pressure  transducers  installed  on  the  upper  surface  of  the  airfoil,  their 
ability  to  capture  shock  location  was  intrinsically  limited  in  resolution.  The 
interferometry  technique  in  their  study  was  utilized  to  study  the  separation  location  for 
different  oscillation  frequencies,  but  was  not  applied  to  the  shock  location. 

Fernie,  Babinsky  and  Bruce  conducted  numerous  investigations  into  the  effects  of 
unsteady  freestream  oscillations  on  the  dynamic  shock  structure  on  a  NACA0012  airfoil 
at  low  angles  of  attack.15'19  When  operated  at  transonic  speeds  and  moderately  high 
frequencies  (of  the  order  of  50  Hz),  they  found  that  there  were  differences  in  the  shock 
strength  and  speed,  depending  on  whether  the  freestream  Mach  number  was  increasing  or 
decreasing.  Evidently,  a  decelerating  freestream  caused  a  stronger  shock  to  move  forward 
on  the  airfoil  surface  more  rapidly  than  an  accelerating  freestream  over  the  same  change 
in  Mach  number.  Fernie  and  Babinsky’s  work  utilized  Kulite  pressure  transducers  and 
high  speed  spark  Schlieren  images  to  study  this  fluid  dynamic  phenomenon.  However 
misalignments  in  the  tunnel  sidewalls  and  or  camera  setup  would  cause  differences 
between  tunnel  runs  at  the  same  operating  conditions.  Therefore  the  pressure  taps  data  at 
every  5%  of  the  chord  (c)  were  utilized  for  quantitative  results  resulting  in  error  bars  of 
+/- 2.5%  x/c. 

At  a  certain  regime  of  steady  Mach  numbers  and  angles  of  attack,  airfoils 

encounter  a  self  induced  phenomenon  known  as  buffeting  where  the  standing  shock  wave 

undergoes  oscillations  in  the  chordwise  direction  of  the  airfoil,  as  seen  in  a  schematic 

diagram  from  Lee,  Figure  5.  Buffeting  has  been  described  by  Crouch  et  al.  as  large 

scale  oscillations  that  can  limit  the  flight  envelope  of  aircraft  due  to  fluid-structure 

interactions  and  are  encountered  by  a  global  instability  in  the  flow  field  that  is  self 
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sustaining.'  The  fluctuations  in  coefficient  of  pressure  ( Cp)  have  been  measured  on  the 

order  of  0.5  due  to  the  fluctuation  in  the  shock  location  over  the  given  pressure  tap  with 

20 

an  oscillation  period  of  ~  14  ms.' 
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Figure  5:  Schematic  diagram  of  self-sustained  shock  oscillations' 

In  order  to  develop  a  full,  accurate  understanding  of  this  dynamic  stall 
phenomenon,  the  problem  must  be  investigated  with  precise  high  speed  instrumentation 
in  conditions  that  adequately  model  the  time-varying  compressibility  effects.  The  purpose 
of  this  investigation  is  to  experimentally  model  unsteady  compressible  flow  and  study  the 
complex  compressibility  effects  on  shock  movement  at  low  reduced  frequencies, 
compressible  Mach  numbers,  and  high  Reynolds  numbers. 

Wind  tunnel  testing  and  computational  fluid  dynamics  (CFD)  have  both  become 
integral  tools  in  the  design,  testing,  and  evaluation  of  airfoils.  Static,  or  time-invariant, 
tests  are  ordinarily  conducted,  which  has  become  commonplace  in  the  aerodynamic 
community.  Static  testing  is  significantly  easier  to  model  than  dynamic  testing  and 
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requires  substantially  less  computational  time  for  CFD.  However  with  the  advancement 
in  computational  power,  high  fidelity  CFD  models  and  wind  tunnel  measurement 
techniques  alike  are  now  possible.  These  methods  can  be  applied  to  perfonn 
comprehensive  unsteady  testing  and  evaluation  of  airfoils.  The  most  productive  approach 
is  one  that  combines  experimental  and  numerical  modeling.  Therefore,  despite  the 
scarcity  of  high-Reynolds  number,  compressible,  dynamic  wind  tunnels,  such  facilities 
are  necessary  to  completely  understand  the  relevant  flow  physics  and  validate  the  CFD 
modeling.  For  these  reasons  The  Ohio  State  University  upgraded  the  6"  x  22"  Transonic 
Wind  Tunnel  to  create  unsteady  compressible  flow  field  as  discussed  in  Chapter  2. 
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Chapter  2:  Freestream  Mach  Number  Oscillation  Mechanism  Design 


Section  2.1:  Freestream  Mach  Number  Oscillation  Introduction 

A  freestream  oscillation  mechanism  is  required  to  create  an  unsteady  flow  field 
for  testing.  There  has  been  some  fundamental  work  in  this  area  of  research  in  the 
incompressible  regime.  Shih  and  Ho  used  a  low  speed  unsteady  water  tunnel  to 
investigate  the  evolving  vorticity  field  about  an  NACA  0012  at  the  incompressible  static 
stall  angle  (a  =  12°)  in  an  unsteady  sinusoidal  freestream  (low."  They  were  able  to  see  an 
increase  in  lift  that  was  co-located  with  a  vortex  being  shed  from  the  airfoil  surface  using 
two  axis  load  cells  and  flow  visualization.  They  utilized  a  vertical  mounted  water  tunnel 
whose  freestream  speed  was  precisely  controlled  by  attaching  a  stepper  motor  to  a  flow 
control  gate  downstream  of  the  test  section.  In  this  manner  different  velocity  profiles 
were  achieved  such  as  saw-tooth  and  sinusoidal  waveforms.  Saxena  performed  similar 
low  speed  oscillating  flow  experiments  on  a  NACA  0012.  This  experimental  setup  at 
Illinois  Institute  of  Technology  oscillated  a  set  of  4  flat  plates  at  the  aft  end  of  the  3-30 
m/s  wind  tunnel  to  create  the  required  flow  conditions.  Favier  et  al.  created  a  more 
unusual  mechanism  by  oscillating  the  airfoil  in  the  streamwise  direction  and  additionally 
in  pitch.9  Other  methods  for  oscillating  the  freestream  flow  include  the  use  of  a  splitter 
plate,  rotating  the  blockage  flow  downstream  or  upstream  of  the  test  section  and 
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oscillating  the  model  relative  to  the  freestream.24  However  the  majority  of  compressible 

13  15  25  26 

flow  facilities  utilize  one  or  more  vanes  rotated  downstream  of  the  test  section.  ’ 

The  6”  x  22”  Transonic  Wind  Tunnel  has  recently  been  modified  to  enable  freestream 
Mach  oscillations  in  order  to  simulate  compressible  dynamic  stall. 


Section  2.2:  Experimental  Facility 

Ohio  State’s  6”  x  22”  Transonic  Wind  Tunnel  was  designed  specifically  for  the 
testing  of  airfoils.  Starting  in  1976  the  6”  x  22”  Transonic  Wind  Tunnel  was  utilized  as 
the  prime  facility  for  The  General  Aviation  Airfoil  Design  and  Analysis  Service.  During 
this  period  of  time  a  plethora  of  airfoils  were  tested  and  developed  at  The  Ohio  State 

27  28 

University  for  the  aeronautical  community,  see  Figure  6.“  ’ 
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High  Pressure  Air 
(42.5  m3  at  15.5MPa) 


Figure  6:  Schematic  diagram  of  the  Aeronautical  &  Astronautical  Research  Laboratory 
6”  x  22”  Transonic  Wind  Tunnel 


The  wind  tunnel  was  designed  as  a  low-turbulence  facility  with  the  primary  goal 
of  matching  freestream  flight  conditions.  To  accomplish  this  objective,  the  settling 
chamber  is  equipped  with  a  perforated  plate,  honeycomb  section  and  then  five  screens  of 
60  mesh  to  lower  the  turbulence  to  acceptable  levels,  less  than  0.5%  under  steady  flow 
conditions.'  A  subsonic  nozzle  with  a  contraction  ratio  of  15:1  provides  excellent  flow 
uniformity  in  the  1.1  m  long,  6”  x  22”  test  section.  The  56  cm  (22  in.)  high  solid 
sidewalls  hold  the  airfoil.  The  15  cm  (6  in.)  spanwise  floor  and  ceiling  walls  are 
perforated  with  3.2  mm-diameter  straight  holes,  yielding  an  effective  porosity  of  6 
percent.  These  isolation  cavities  or  plenums  are  open  to  the  flow  only  at  the 
downstream  end  of  the  test  section  and  aid  in  producing  a  high-quality  flow  in  the  test 
section  by  reducing  Mach  wave  reflections  in  transonic  flow. 
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The  tunnel  is  fed  by  two  21.25  m  tanks  pressurized  up  to  15.5  MPa  with 
conventional  air  dryers  used  to  maintain  gas  purity  with  a  dew  point  of  -40°C.  The  high 
pressure  air  is  controlled  by  two  valves.  The  first  is  a  control  valve  which  sets  the  total 
pressure  and  Reynolds  number.  The  second  is  a  fast-acting  valve  used  to  start  and  stop 
the  flow.  Mach  number  is  controlled  independently  of  Reynolds  number  by  adjusting  the 
throat  area  downstream  of  the  test  section.  The  Mach  number  is  uniquely  set  by  the  throat 
area,  independent  of  stagnation  pressure,  as  long  as  choked  flow  is  maintained.  The 
minimum  Reynolds  number,  at  a  given  Mach  number,  is  set  by  the  stagnation  pressure 
needed  to  choke  the  flow  while  the  maximum  is  set  by  a  structural  pressure  limitation  of 
the  tunnel.  These  two  parameters  (Re  and  M)  can  be  independently  varied  over  a 
considerable  range,  as  shown  in  Figure  7.  The  run  duration  is  typically  10  seconds  with 
15  minutes  between  runs,  however  run  durations  up  to  45  seconds  are  possible  in  this 
blowdown  facility. 
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Figure  7:  Range  of  test  section  operating  conditions  for  the  6”  x  22”  Transonic  Wind 
Tunnel  at  280K  and  a  maximum  operating  pressure  of  345kPa 


The  Mach  number  in  the  test  section  is  detennined  by  the  Mach  -  area  equation 

•  •  9  #  #  # 
seen  in  Equation  1.  Where  A  is  the  cross  sectional  area  at  the  choke  point  and  A,  is  the 

cross  sectional  area  in  the  test  section.  Similar  to  other  wind  tunnel  facilities  as  discussed 

in  Section  2.1,  the  freestream  Mach  oscillation  was  created  by  varying  the  throat  area 

downstream  of  the  test  section  periodically. 
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Equation  1 :  Mach  -  area  relation 
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Section  2.3:  Sinusoidal  Vane  Design 

A  rectangular  vane  set  was  initially  implemented  in  the  wind  tunnel. ~  Even 
though  the  resulting  Mach  number  variation  was  expected  to  contain  higher  hannonics, 
this  geometry  was  selected  for  ease  of  manufacturing  and  implementation.  The  results 
from  the  rectangular  vanes  were  analyzed  to  aid  in  the  design  of  new  sinusoidal  Mach 
number  oscillation  vanes.  Each  of  the  four  rectangular  vanes  had  a  length  of  6.35  cm  and 
a  width  of  3.81  cm.  The  rectangular  vane  set  was  implemented  in  the  tunnel  down-stream 
of  the  test  section  as  seen  in  Figure  6  labeled  “choke  vanes.”  The  tunnel  was  then  run  at  a 
modest  Mach  oscillation  frequency  (/)  of  2.5  Hz.  The  total  (/;„)  and  static  pressures  (/;,) 
were  measured  in  the  test  section  in  order  to  characterize  the  resulting  Mach  oscillation. 
Figure  8  shows  the  measured  and  theoretical  test  section  Mach  number  as  derived  from 
Equation  1,  using  the  Mach  -  area  relation.  Similarly,  the  actual  throat  area  (A*)  is 
calculated  by  solving  Equation  1  with  the  experimentally-measured  test  section  Mach 
number  and  test  section  area. 
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Figure  8:  Experimental  (f=  2.5  Hz)  and  theoretical,  A*  and  Mach  number  for  one  period 
of  the  flow  oscillation  with  rectangular  choked  vanes 


Figure  8  illustrates  that  the  measured  Mach  number  follows  the  general  trend 
anticipated  from  the  Mach  -  area  relation  (Equation  1),  but  with  rounding  of  sharp 
comers  in  the  waveform.  This  attenuation  of  higher  frequencies  is  due  to  the  low-pass 
filter  characteristics  of  the  wind  tunnel,  which  are  discussed  in  more  detail  in  Section  3.2. 
The  difference  between  the  measured  and  theoretical  Mach  number  results  is  largest 
when  the  axes  of  the  rectangular  vanes  are  parallel  with  the  wind  tunnel  axis  (0  and  90 
degrees).  As  the  boundary  layer  increases  along  the  walls  of  the  tunnel  the  displacement 
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thickness  is  found  in  order  to  calculate  the  effective  cross  sectional  area  at  the  test  section 


and  throat  location.  The  throat  area  has  a  larger  displacement  thickness  due  to  the  throat 
location  being  downstream  of  the  test  section,  the  effective  area  ratio,  A  /At,  is  reduced 
giving  a  proportionally  lower  Mach  number  as  seen  in  Equation  1.  The  actual,  or 
experimentally  measured,  Mach  number  is  lower  from  the  theoretical  values  due  to  these 
viscous  effects,  as  the  boundary  layer  along  the  tunnel  walls  grow  from  the  test  section  to 
the  throat. 

Due  to  difficult  to  account  for  boundary  layer  and  viscous  effects  in  the  tunnel 

•  •  •  •  •  •  ♦ 
and  on  the  Mach  oscillation  vanes,  an  empirical  relationship  between  M  and  A  was 

defined.  The  viscous  effects  that  cause  deviation  between  theory  and  experiment  are 

difficult  to  measure;  thus,  an  empirical  relationship  between  Mach  number  and  vane  area 

is  used  for  subsequent  vane  designs.  Figure  9  shows  the  rectangular  choke  vane  static 

results  obtained  with  the  vanes  rigidly  fixed  in  the  0  and  90  degree  orientations.  Results 

from  the  dynamic  oscillation  (Figure  8)  are  also  indicated  in  Figure  9,  after  low-pass 

filtering  to  remove  the  higher  hannonics.  The  static  results  correlate  well  to  the  dynamic 

/=  2.5  Hz  data.  A  second  order  polynomial  fit  (see  Equation  2)  was  applied  to  the 

experimental  data  of  Figure  9  to  detennine  the  empirical  relationship  between  the  test 

section  Mach  number  and  throat  area.  In  addition  to  the  data  points  plotted,  M  =  0  and  1 

were  also  constrained  as  full  and  zero  blockage  area,  respectively. 
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Figure  9:  Empirical  relationship  between  the  test  section  Mach  number  and  throat  area 
for  use  in  sinusoidal  choke  vane  design 


Equation  2:  Empirical  relationship  between  Mach  number  and  throat  area. 
21*  =  -11.35  M2  +  41.65  M  +  0.4907 


This  empirical  relationship  was  then  used  to  design  an  elliptical  vane  set  that 
would  produce  a  sinusoidal  oscillation  in  Mach  number.  The  vanes  are  designed  with  a 
symmetric  profile  in  order  to  avoid  eccentric  loading;  therefore,  two  periods  of  the  Mach 
wavefonn  are  created  for  each  complete  rotation  of  the  vanes.  The  vane  radius  was 
designed  using  the  empirical  relationship  between  the  test  section  Mach  number  and  the 
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throat  area  defined  in  Equation  2  and  the  following  equation  for  the  desired  Mach 
number: 

Equation  3:  Ideal  Mach  number  oscillation. 

M(t)=  0.105  *sin(t)+ 0.605 

The  resulting  choke  vane  design  for  this  Mach  wave  is  shown  in  Figure  10.  The 
resulting  design  turned  out  to  be  much  too  aggressive  and  vanes  ended  up  peanut  shaped 
because  the  amplitude  of  the  sine  wave  was  too  steep.  Obviously  the  minimum  area  is  not 
fluid  dynamically  possible  to  obtain  using  this  area  ratio  Mach  curve.  The  mean  Mach 
number  and  amplitude  were  reduced  in  order  to  more  closely  resemble  that  of  the 
rectangular  vanes.  This  allowed  for  a  feasible  vane  shape  using  the  A*  -  Mach  curve  fit, 
as  well  as  higher  confidence  in  the  curve  fit  due  to  the  similar  Mach  number  range  of  the 
elliptical  and  rectangular  vane  set.  Figure  9  shows  the  Mach  number  and  A*  for  one 
rotation  of  the  designed  vanes. 
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Length  (cm) 

Figure  10:  Vane  profile  for  a  Mach  Number  =  0.105  *  sin(t)  +  0.605 


The  Mach  oscillation  was  redesigned  based  on  the  new  Mach  number  oscillation  seen  in 
Equation  4.  The  A*  is  again  found  from  inputting  this  Mach  number  equation  into  the  A* 
-  Mach  polynomial  curve  lit.  The  necessary  blockage  area  to  obtain  A*  is  divided 
between  the  four  vanes  and  wrapped  around  the  axis  of  rotation.  This  results  in  the 
constructed  and  tested  vane  geometry  of  Figure  11.  The  set  of  vanes  were  manufactured 
to  a  tolerance  of  0. 127  mm  and  installed  in  the  wind  tunnel  for  the  current  investigation. 
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Equation  4:  Designed  Mach  number  oscillation. 
M{t)  =  0.063  *  sin(t)  +  0.45 


Figure  1 1 :  Choke  vane  profile  constructed  to  produce  sinusoidal  Mach  number  oscillation 


Section  2.4:  Freestream  Oscillation  Mechanism 


The  transonic  wind  tunnel  has  been  modified  to  oscillate  the  flow  using  a 
variable-area  throat  downstream  of  the  test  section.  The  throat  location  can  be  seen  in 
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Figure  6,  labeled  as  the  choke  vanes.  The  choke  vanes  are  connected  to  a  5-HP  variable 
frequency  A/C  motor  through  a  drive  mechanism  (Figure  12).  This  drive  mechanism 
consists  of  a  sprocket  mounted  to  each  vane  with  a  chain  linking  them  to  the  motor.  The 
choke  vanes  are  counter  rotating  to  minimize  harsh  unsteady  vertical  loading  on  the 
tunnel.  An  optical  encoder  is  mounted  to  one  vane  shaft  in  order  to  record  the  rotational 
frequency  of  the  choke  vanes.  This  includes  a  500  per  revolution  and  a  once  per 
revolution  signal  which  is  necessary  to  calculate  the  phase,  compare  waveforms,  and 
produce  a  triggering  mechanism  for  the  particle  image  velocimetry,  pressure-  and 
temperature-sensitive  paint  experimental  setups.  In  the  current  geared  configuration,  the 
A/C  motor  can  produce  flow  oscillatory  frequencies  up  to  21  Hz.  This  oscillatory 
frequency  is  important  in  its  relation  to  the  reduced  frequency  (k),  as  seen  in  Equation  5. 
For /=  21  Hz,  M=  0.52  and  c  =  0.127  m,  k  =  0.05,  which  is  traditionally  defined  as  the 

3 1 

end  of  quasi  steady  and  beginning  of  unsteady  flow. 

Equation  5:  Reduced  frequency  equation. 

k  =  £ 

2V 
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Figure  12:  Mach  oscillation  drive  train  assembly 
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Chapter  3:  Unsteady  Mach  Oscillation  Characterization 


Section  3.1:  Wind  Tunnel  Characterization  Experimental  Setup 

A  TSI  model  1231-60W  hot  film  was  the  main  measurement  instrument  used  for 
the  initial  characterization  tests.  The  hot  film  was  used  to  investigate  the  effects  of 
changing  Mach  oscillation  frequencies  on  the  wavefonn  shape  and  phase  shift.  The  hot 
film  is  surface  mounted  on  a  conical  probe  which  is  centered  in  the  test  section.  A  thin 
platinum  film  strip  of  0.127  mm  width  is  wrapped  around  the  conical  probe  tip,  and  the 
tip  is  aligned  to  point  upstream.  The  hot-film  signal  was  sampled  at  2  kHz  and  low-pass 
filtered  at  1  kHz  to  attenuate  noise  and  prevent  aliasing. 

Measurements  were  made  of  the  magnitude  and  phase  relationship  of  the 
oscillating  flow  through  the  wind  tunnel  at  unsteady  Mach  numbers  ranging  between  0.44 
<  M  <  0.65  and  at  various  Mach  number  oscillation  frequencies  (f)  up  to  21  Hz.  During  a 
30-second  run,  the  total  temperature,  T0,  typically  decreases  15  K,  corresponding  to  a  7.5 
m/s  decrease  in  the  speed  of  sound  (a).  Due  to  large  tubing  effects  on  the  pressure  tap 
lines  and  low  sampling  frequencies  of  the  pressure  transducers  used  in  this  initial  setup, 
the  maximum  oscillation  frequency  used  was  2.1  Hz.  These  low  frequency  Mach 
oscillation  data  from  pressure  transducers  were  utilized  to  perform  an  in  situ  calibration 
of  the  hot  film  by  throttling  back  the  oscillation  frequency  to  2.1  Hz  at  the  beginning  and 
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end  of  each  run.  The  combined  uncertainty  in  measured  Mach  number  is  ±  0.005  at  a 
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nominal  Mach  number  of  0.51.  " 


Section  3.2:  Unsteady  Tunnel  Characterization  Results  and  Discussion 

The  sinusoidal  vanes  for  the  geometry  seen  in  Figure  1 1  were  manufactured  to  a 
geometrical  tolerance  of  0.127  mm  and  installed  in  the  wind  tunnel.  A  series  of  wind 
tunnel  runs  were  completed  to  analyze  their  performance.  During  characterization  the 
vanes  were  oscillated  at  frequencies  ranging  from  2.5  to  21  Hz.  The  sinusoidal  choke 
vanes  were  found  to  produce  a  periodic  Mach  oscillation  from  0.44  to  0.64  with  no  airfoil 
in  the  test  section.  The  maximum  Mach  number  was  shifted  up  from  the  designed 
Equation  4  value  by  AM  =  0.127  and  the  minimum  by  AM  =  0.053.  The  sinusoidal  data  is 
seen  to  be  closer  to  that  of  the  A*/A  (M)  equation.  Equation  1,  in  Figure  13.  The  Mach 
oscillation  waveform  is  seen  (Figure  14)  to  be  bottom  heavy  due  to  the  difference 
between  the  instantaneous  angular  radius  and  the  maximum  thickness  for  the  given 
angular  position.  The  Mach  number  also  had  a  bulk  shift  of  the  mean  Mach  number  from 
this  instantaneous  angular  radius  mean  M=0.54  to  the  experimentally  measured  M=0.515, 
as  seen  by  the  horizontal  lines  in  Figure  14.  Therefore  the  1-D  analytical  Mach  number 
accounts  for  the  majority  of  these  discrepancies  in  instantaneous  cross  sectional  area  by 
including  area  changes  due  to  boundary  layer  thickness  at  the  test  section  and  choke  area, 
as  discussed  by  Gompertz  et  ah.  The  experimentally  measured  Mach  number  is  still 
lower  than  the  1-D  analytical  solution,  effectively  stating  that  the  throat  area  of  the  tunnel 
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is  smaller  than  predicted.  It  is  postulated  that  this  added  throat  blockage  area  is  due  to 
unsteady  wakes  of  the  rotating  vanes  or  non-uniform  boundary  layer  growth  on  the  vanes. 
Figure  14  shows  that  the  rectangular  vanes  experienced  additional  blockage  possible  due 
to  larger  wakes  shedding  from  the  discontinuous  rectangular  vanes.  The  experimentally 
measured  Mach  oscillation  wavefonn  is  also  slightly  bottom  heavy,  with  the /=  2.5  Hz 
mean  Mach  number  4.5%  lower  than  the  peak  to  peak  median.  Figure  14  shows  both  the 
1-D  analytical  designed  Mach  oscillation  as  well  as  the  Mach  wave  achieved  in  the  test 
section  for  the  sinusoidal  vane  set  at  a  frequency  of  2.5  Hz. 


Figure  13:  Throat  Area  vs  Mach  number  for  sinusoidal,  rectangular  static  and  dynamic 
data  as  well  as  Equation  1  and  the  designed  curve  fit 
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Figure  14:  Designed  test  section  Mach  number  oscillation  compared  to  the  test  section 
Mach  number, /=  2.5  Hz33 

Equation  6  is  the  1-D  wave  propagation  speed.  This  is  significant  for  this  tunnel 
because  it  is  operated  in  the  compressible  Mach  regime  and  the  test  section  is  1.75  m 
upstream  from  the  choke  vanes.  As  the  Mach  number  increases,  the  flow  requires  a 
longer  amount  of  time  to  propagate  up  the  tunnel.  The  rectangular  vane  data  seen  in 
Figure  9  shows  this  hysteresis  effect  as  the  data  does  not  collapse  on  itself  for  the  upward 
and  downward  slope  of  Mach  number.  Therefore,  by  the  time  the  Mach  wave  reaches  the 
test  section,  the  wave  will  be  shifted  with  the  higher  Mach  number  flow  lagging  the  lower 
speed  flow. 

Equation  6:  Wave  propagation  speed. 

C  =  a(l  -  M) 

A  phase  shift  and  wave  steepening  effect  is  observed  as  the  period  of  the  wave  decreases 
while  the  propagation  speed  remains  the  same.  Figure  15  shows  how  the  waveform 
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changes  with  different  frequencies,  the  Mach  oscillation  waves  are  shifted  in  phase  to 
align  at  the  maximum  Mach  number.  It  was  also  a  prerequisite  to  normalize  the  Mach 
number  due  to  the  amplitude  attenuation  at  the  higher  frequencies.  The  wavefonn  shape 
is  seen  to  change  from  a  bottom  heavy  shape  seen  with  the/=2.5  Hz  oscillation  to  a  more 
symmetrical  sine  wave  in  the  /  =16  Hz  oscillation.  The  wave  steepening  from  the 
propagation  speed  can  be  seen  in  Figure  15,  with  the  minimum  normalized  Mach  number 
moving  to  an  earlier  position  in  the  period  as  the  frequency  is  increased. 


Figure  15:  One  period  of  normalized  hot  film  data  plotted  at  different  frequencies 
illustrating  wave  steepening 


The  Mach  number  oscillation  was  further  tested  at  frequencies  up  to  21  Hz.  The 
hot  wire  study  was  able  to  characterize  the  Mach  number  oscillation  amplitude  and  phase 
delay  for  this  range  of  frequencies.  The  amplitude  and  phase  shift  bode  plots  of  Mach 
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number  oscillation  are  shown  in  Figure  16.  The  oscillation  amplitude  is  attenuated 
significantly  at  high  oscillation  frequencies.  Bode  plots  indicate  that  the  -3  dB  roll-off 
point  is  at  approximately  8  Hz  for  an  empty  test  section.  The  phase  shift  is  a  linear 
function  of  Mach  number  oscillation  frequency,  for  a  given  mean  Mach  number,  which 
indicates  that  the  wavespeed  is  approximately  constant.  Figure  17  illustrates  the  phase 
shift  and  amplitude  attenuation  for  four  representative  Mach  oscillation  waveforms  at 
different  frequencies. 


Mach  Number  Oscillation  Frequency,  f(Hz) 


Figure  16:  (a)  Amplitude  and  (b)  Phase  Number  vs.  Mach  Number  Oscillation  Frequency 

if)  33 
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Mach  Number, 


Figure  17:  Representative  test  section  Mach  number  waveforms  at  different  oscillation 

-5-5 

frequencies 


Chapter  4:  Pressure  Sensitive  Paint 

Section  4.1:  Pressure  Sensitive  Paint  Introduction 

Pressure-sensitive  paint  (PSP)  is  an  experimental  technique  based  on  the 
quenching  mechanisms  of  luminescence  molecules.  Compared  with  conventional 
techniques  like  pressure  taps,  it  offers  a  unique  capability  for  high  frequency,  non- 
contact,  full-field  measurements  of  surface  pressure  on  a  complex  aerodynamic  model 
with  higher  spatial  resolution  at  a  lower  cost.  PSP  has  to  be  excited  by  a  light  source  to 
trigger  the  luminescence  response,  which  is  the  paint  down-converting  the  excitation 
illumination  to  a  longer  wavelength.  The  intensity  of  this  emitted  light  is  inversely 
proportional  to  the  amount  of  oxygen  in  the  air  and  is  also  a  function  of  the  surface 
temperature  of  the  model;  as  the  luminescent  intensity  (I)  of  the  phosphor  is  determined 
by  the  thermal  and  oxygen  quenching  mechanisms  of  the  PSP.  Temperature-sensitive 
paint  (TSP)  is  insensitive  to  oxygen  concentration,  so  its  luminescent  intensity  is  only  a 
function  of  the  model  temperature. 
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Section  4.2:  Lifetime-  vs.  Intensity-Based  PSP  Techniques 


Lifetime-  and  intensity-method  are  the  two  major  methods  of  recording  the 
luminescent  emission  response  from  the  PSP.  The  first  method  is  lifetime-based,  where 
the  PSP  is  exposed  to  a  pulsed  excitation  light.  A  double  shutter  scientific  camera  then 
images  the  luminescent  intensity  during  the  excitation  illumination  and  the  corresponding 
lifetime  decay  (x)  of  the  paint,  which  is  shown  as  Equation  7  and  illustrated  in  Figure 
18.34'36  The  lifetime -method  for  PSP  is  based  on  the  fact  that  the  luminescence  decay 
lifetime  is  inversely  proportional  to  the  local  oxygen  concentration,  which  is  directly 
related  to  the  air  pressure  according  to  Henry’s  law.  This  relation  is  shown  as  the  Stem- 
Volmer  equation,  Equation  8,  where  xref  and  pref  are  the  luminescent  lifetime  and  air 
pressure  at  a  reference  condition,  respectively.  The  Stem-Volmer  coefficients  A  and  B, 
which  are  temperature-dependent  due  to  thennal-quenching,  are  experimentally 
detennined  by  calibration. 

Equation  7:  The  luminescent  emission  response. 

t 

l(t)-Ae  r 

Equation  8:  Stern- Volmer  equation. 

^-  =  A(T)  +  B(T)-^- 

r  Pref 
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Figure  18:  Representation  of  two-gate  lifetime  PSP  measurement  technique37 


Lifetime-based  PSP  measurement  approaches  include  phase-sensitive  detection 
and  multi-gate  integration  techniques.  A  typical  schematic  of  the  two-gate  lifetime 
approach  is  shown  in  Figure  18.  In  the  lifetime-based  technique,  paint  excitation  is 
provided  by  a  modulated  light  source  such  as  a  pulsed  laser  or  LED  flash  lamps. 
Generally,  paint  emission  intensities  are  acquired  and  integrated  over  two  distinct  gates, 
where  the  first  gate  occurs  during  the  illumination  pulse  and  the  second  is  centered  on  the 
decay  of  the  fluorescence.  In  this  application,  however,  two  nearly  contiguous  gates  are 
used,  based  on  the  double  shutter  mode  of  a  charge-coupled  device  (CCD)  camera.  The 
first  gate  (Gl)  captures  the  signal  in  the  illumination  pulse,  while  the  second  gate  (G2) 
captures  the  lifetime  decay.  The  integrated  signal  of  the  first  gate  is  relatively  insensitive 


to  pressure,  thus  serving  as  a  reference,  while  the  second  gate  signal  is  sensitive  to 
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pressure  changes  according  to  the  Stern-Volmer  equation.  The  ratio  of  G1  and  G2  is  able 
to  provide  surface  pressure  information  after  calibration.  In  practice,  the  timing  of  these 
gates  to  optimize  the  pressure  sensitivity  of  the  lifetime  PSP  method  and  signal-to-noise 
ratio  (SNR)  is  found  at  equal  intensities  of  G1  and  G2.  The  major  advantage  of  the 
lifetime  approach  is  the  minimization  of  illumination  errors  by  self  referencing  G2  with 
G1  for  both  the  wind-on  and  wind-off  images. 

Theoretically,  the  lifetime  of  the  paint  at  uniform  pressure  and  temperature  should 
be  constant  at  each  point  on  the  surface.  However,  it  has  been  shown  by  Goss  et  al.  that 
this  is  not  the  case.  ’  This  spatial  variation  may  be  due  to  non-uniform  illumination 
field  or  variations  in  luminophore  concentration.  While  the  source  of  this  noise  remains 
unclear,  it  has  been  demonstrated  that  a  wind-off  ratio  can  be  used  to  minimize  this  noise 
in  both  the  lifetime  and  intensity  methods.  In  addition  to  illumination  errors  the 
temperature  sensitivity  of  PSP  is  a  source  of  uncertainty  (see  Section  6.5). 

When  the  lifetime-based  method  is  applied  using  a  single-shot  technique,  it  is 
capable  of  eliminating  the  excitation  variations,  model  movements  and  model 
deformations  due  to  the  technique  being  self-referenced.  In  this  method,  a  high-power 
laser  is  selected  as  the  illumination  source,  for  the  purpose  of  improving  signal  to  noise 
ratio  during  data  acquisition.  Due  to  the  pulse-to-pulse  non-uniformity  of  laser 
illumination  it  is  always  necessary  to  utilize  the  lifetime  mode,  where  both  the  pressure 
and  reference  images  are  obtained  on  the  same  decay  curve  from  a  single  laser  pulse. 
This  single-shot  lifetime  PSP  method  was  applied  by  Kumar  for  surface  pressure 
measurements  on  rotating  surfaces  with  unsteady  pressure  fluctuations.40  It  is  also 


suitable  for  unsteady  applications  with  large  pressure  gradients  and  model  movements, 
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such  as  oscillating  airfoils  in  compressible  flow.  Additionally  the  single-shot  lifetime- 
method  was  utilized  in  the  6”  x  22”  Transonic  Wind  Tunnel  for  the  initial  Mach 
oscillation  PSP  study  by  Jensen  et  ah,  the  schematic  diagram  of  the  setup  can  be  seen  in 
Figure  19. 41  This  study  utilized  a  fast-response  porous  polymer  paint  which  was  first 
developed  by  Scroggin  et  ah. 42  This  paint  is  a  polymer/ceramic  based  PSP  which  was 
designed  specifically  to  improve  the  diffusivity  of  the  binder  material,  thus  allowing 
faster  response  times.  Its  kinetic  characteristics  were  studied  in  detail  by  Gregory  et  ah;  it 
has  a  flat  response  time  greater  than  5  kHz,  which  is  sufficient  to  study  the  flow  field  on 
airfoils  with  oscillating  freestream  frequencies  less  than  21  Hz.43,  44  PtTFPP,  a  typical 
platinum  porphyrin-based  luminophore,  was  used  as  the  PSP  sensor  in  the  initial 
unsteady  PSP  study  by  Jensen  et  al.  and  as  the  pressure  channel  of  the  bi-luminophore 
paint  during  the  later  NACA  0012  steady  and  unsteady  experiemnts.45 
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Figure  19:  Schematic  diagram  of  the  single-shot  lifetime-method  PSP  experimental  setup 


The  other  major  PSP  method  is  called  the  intensity-based  method,  wherein  the 
illumination  source  is  continuously  on  and  the  camera  records  the  intensity  of  paint.  This 
requires  the  illumination  source  to  be  very  consistent.  It  has  the  advantage  of  being  a 
much  simpler  technique  and  does  not  require  a  double  shutter  mode  camera.  The  phase 
locked  triggering  of  the  camera  shutter  is  still  onerous  yet  simpler  when  compared  to  that 
of  the  lifetime -method,  due  to  the  lifetime-based  method  requiring  specific  timing 
between  the  illumination  pulse  and  both  camera  shutters.  The  drawback  of  the  intensity- 
based  method  is  that  continuous  and  consistent  illumination  provided  by  LED’s  is  not  as 
powerful  as  with  a  pulsed  laser.  In  the  past  this  has  caused  the  need  for  phase  averaging 
to  boost  the  SNR  of  the  paint  to  acceptable  levels.  In  recent  years,  however,  illumination 
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sources  have  become  powerful  enough  to  acquire  all  of  the  data  needed  in  a  single  shot, 
similar  to  the  method  described  with  the  lifetime  method.  Initial  experiments  on  a  NACA 
0021  utilized  the  lifetime  method  with  400  mJ  per  pulse  of  energy.  However,  the 
Nd:YAG-532  mn  laser  as  an  illumination  source  is  not  compatible  with  the  bi- 
luminophore  paint  because  it  overlaps  with  the  emission  spectrum  of  the  temperature 
channel.  Therefore,  a  set  of  three  LEDs  were  utilized  to  provide  18  W  of  continuous  405- 
nm-wavelength  illumination  on  the  airfoil  surface.  This  advance  in  LED  illumination 
power  effectively  enables  the  pressure  and  temperature  sensitive  paint  to  make  single¬ 
shot  measurements.  The  single-shot  intensity-method  was  chosen  for  experimental  setup 
for  its  ability  to  instantaneously  acquire  the  unsteady  temperature  and  pressure  map. 


Section  4.3:  Temperature  Compensating  Pressure  Sensitive  Paint 

The  single-shot  pressure-sensitive  paint  technique  is  new  compared  to  the  more 
common  phase-averaging  PSP  technique.40  The  advantage  of  the  single-shot  PSP  for 
unsteady  measurements  is  the  ability  to  capture  all  the  data  within  a  few  microseconds. 
Transient  effects  that  would  normally  be  averaged  out  under  the  phase-averaging 
technique  are  captured  with  the  single-shot  technique.  Jensen  et  al.  tested  the  validity  of 
this  method  in  a  changing  temperature  environment.41  A  NACA  0021  airfoil  was  tested  in 
the  6”  x  22”  Transonic  Wind  Tunnel  in  oscillatory  flow  to  evaluate  this  PSP  technique 
for  making  unsteady  pressure  measurements  and  to  begin  study  of  the  transient  shock 
characteristics  in  quasi-steady  flow.  The  airfoil  was  tested  at  a  Mach  oscillation 
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frequency  of  2.45  Hz.  An  example  coefficient  of  pressure  (Cp)  profile  plot  and  upper 
surface  Cp  contour  map  can  be  seen  in  Figure  20.  Due  to  the  wind  tunnel  being  a 
blowdown  facility,  the  stagnation  temperature  drops  throughout  the  run,  typically  15  K 
during  a  25  second  run.  These  temperature  effects  on  the  PSP  can  be  seen  as  a  false  rise 
in  the  Cp  of  the  airfoil  at  the  trailing  edge  and  the  spanwise  non-unifonnity  in  the  upper 
surface  Cp  map,  as  shown  by  a  theoretical  black  dashed  line  in  Figure  20.  The  false  rise  is 
due  to  a  temperature  change  from  the  wind-on  to  the  wind-off  images  and  is  intrinsic  to 
the  PSP’s  temperature  sensitivity.  An  in  situ  calibration  is  able  to  compensate  for 
variations  in  bulk  temperature  shifts;  however  spanwise  and  chordwise  variations  in 
temperature  will  still  be  prevalent  sources  of  error  in  the  data.  The  chordwise  temperature 
variation  can  be  seen  most  dramatically  at  the  aft  end  of  the  airfoil,  where  it  is  the 
thinnest  and  has  the  largest  gradients  in  temperature.  Additionally  the  spanwise  variations 
are  due  to  the  model  conducting  heat  through  the  tunnel  sidewalls  due  to  different 
endplate  shapes  on  the  airfoil  model. 
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Figure  20:  Single-luminophore  PSP  derived  Cp  plot  (left)  and  upper  surface  Cp  map 
(right)  for  a  NACA  0021  under  oscillatory  flow  at  f=2. 5  Hz,  a=10°,  Re  =4.1  million,  and 
M=0.61 


In  this  lifetime -based  test  the  PSP  was  calibrated  via  an  in  situ  method  using  the 
pressure  taps  at  x/c  of  15  and  50%  on  the  upper  surface.  Only  two  pressure  taps  are 
needed  because  the  pressure  calibration  is  assumed  to  be  linear  over  this  small  pressure 
range.  Furthennore,  the  accuracy  of  the  calibration  over  the  first  half  of  the  airfoil  is 
illustrated  by  how  well  the  third  and  fourth  pressure  taps  agree  with  the  calibrated  PSP,  as 
seen  by  the  red  circles  in  Figure  20  (left).  This  in  situ  calibration  is  necessary  to  account 
for  the  temperature  drop  throughout  the  operation  of  the  blowdown  wind  tunnel. 
However,  this  technique  only  provides  a  bulk  correction  for  temperature  effects  across 
the  solid  aluminum  airfoil.  This  effect  is  shown  in  Figure  20,  where  the  slope  of  Cp 
changes  past  x/c  of  70%.  Even  with  the  poor  calibration  at  the  trailing  edge,  this  figure 
illustrates  the  detail  PSP  can  offer  in  spatial  resolution  and  the  detail  in  which  a  relatively 
weak  standing  shock  wave  can  be  resolved. 
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Due  to  these  large  variations  in  temperature  throughout  the  tunnel  run  a 
temperature  correction  method  is  needed  in  order  to  make  accurate  PSP  measurements, 
which  is  discussed  further  in  Section  6.4.  The  method  used  a  bi-luminophore  pressure- 
and  temperature-sensitive  paint.  A  typical  pressure  sensor  and  temperature  sensor  were 
combined  into  one  paint  fonnulation  and  were  separated  spectrally  with  regards  to  their 
emission  signals.  This  bi-luminophore  method  allows  for  the  addition  of  a  reference 
channel  that  is  only  sensitive  to  temperature  to  account  for  the  pressure  channel’s 
temperature-sensitivity.  Thus,  an  accurate  PSP  measurement  is  obtained  while  also 
acquiring  the  temperature  map  for  the  model  if  needed.  Various  other  authors  have 
attempted  similar  bi-luminophore  techniques  as  summarized  in  Peng  et  al..46’  47 
Alternative  temperature  correction  techniques  have  used  infra-red  cameras  and  bifurcated 
airfoils  painted  half  with  each  TSP  and  PSP. 


Section  4.4:  Two  Camera  vs.  Color  Camera  PSP  and  TSP  Methods 

Two  different  camera  methods  were  evaluated  for  this  temperature  compensated 
PSP  measurement.  A  standard  color  camera  with  a  built-in  red  green  blue  (RGB)  Bayer 
filter  over  the  individual  pixels  was  initially  employed.46,  47  By  only  having  one  camera 
this  technique  was  very  easy  to  implement  optically.  It  was  found  that  the  Bayer  filter  on 
the  camera,  an  Alta  U2000C,  had  significant  filter  overlap  between  the  color  channels, 
which  added  pressure  sensitivity  to  the  temperature  channel.  It  was  also  challenging  to 
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create  a  bi-luminophore  paint  formulation  that  had  similar  intensities  for  both  the 
pressure  and  temperature  channels  at  all  operating  conditions. 

The  two  camera  method  was  therefore  investigated,  compromising  the  ease  of 
optical  setup  for  that  of  flexibility  to  optimize  the  intensity  of  each  camera  individually 
by  means  of  pixel  binning  and  aperture  settings.  Additionally  individual  filters  could  be 
chosen  for  each  black  and  white  PCO-1600  camera  in  order  to  minimize  the  bleed 
through  of  each  color  channel  and  give  better  optical  spectral  separation.  For  these 
reasons  the  two-camera  method  was  chosen  for  this  experiment.  A  direct  comparison  of 
the  unfiltered  contours  of  Cp  measured  with  the  two  methods  can  be  seen  in  Figure  2 1 . 
The  two-camera  method  is  seen  to  have  less  noise  and  thus  a  better  signal  to  noise  ratio. 


Color  Camera  Technique  Two  Camera  Technique 


Figure  2 1 :  Comparison  of  the  raw  temperature-corrected  contours  of  Cp  for  the  Color 
Camera,  Left,  and  the  Two  Camera  technique,  right,  on  a  NACA  0012  at /=  9.5  Hz, 
a=9°,  Re=3.9  million,  and  M=0.57,  flow  from  top  to  bottom 


The  two  methods  were  compared  quantitatively  as  follows.  The  intensity  ratio  of 
the  green  channel  was  recorded  for  both  the  two-camera  and  color-camera  methods  for  a 
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matrix  of  pressures  and  temperatures.  The  calibration  was  performed  in  a  specially  built 
calibration  chamber  designed  to  hold  pressure  and  temperature  in  a  range  from  0  to  241 
KPa  and  -8  to  50°  C,  respectively.  An  Ocean  Optics  USB4000  spectrometer  was  used  to 
illustrate  TSP’s  insensitivity  to  pressure  as  well  as  the  spectral  separation  between  the 
two  channels,  as  seen  in  Figure  22  and  Figure  23,  from  Peng  et  al..47  The  color  and  two- 
camera  filters  were  overlaid  onto  the  spectroscopy  data.  It  can  be  seen  that  the  green  filter 
of  the  color  camera  also  collects  some  of  the  pressure  data  at  700  to  750  nm  wavelengths; 
this  gives  false  pressure  sensitivity  to  the  TSP  measurement.  Figure  24  and  Figure  25 
show  the  temperature  vs.  green  intensity  for  the  color-  and  two-camera  methods, 
respectively,  with  a  reference  temperature  of  20°C.  A  second  order  polynomial  was  fit  to 
the  presented  data  and  characterized  by  the  coefficient  of  determination,  R  .  In  the  color 
camera  data  it  is  apparent  that  the  temperature  channel  is  acquiring  some  of  the 
luminescence  signal  from  the  PSP  thus  showing  some  pressure  sensitivity  due  to  the  poor 
quality  of  the  Bayer  filter.  Surface  fits  were  performed  for  the  pressure  channel  with  their 
respective  R~  values  presented  in  Table  1. 
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A)  Constant  Temperature 


B)  Constant  Pressure 
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Figure  22:  Emission  spectra  for  the  color  camera  technique  at  A)  different  pressures  and 
constant  temperature  and  B)  different  temperatures  and  constant  pressures,  with  the 
respective  green  and  red  band-pass  optical  filters  shown47 


A)  Constant  Temperature 


B)  Constant  Pressure 
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Figure  23:  Emission  spectra  for  the  two-camera  technique  at  A)  different  pressures  and 
constant  temperature  and  B)  different  temperatures  and  constant  pressures,  with  the 
respective  green  and  red  band-pass  optical  filters  shown47 
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Table  1 :  R"  values  for  the  pressure  and  temperature  curve  fits  for  the  color  and  two 
camera  methods. 


Two  Camera 

Color  Camera 

Temperature  R2 

0.998 

0.978 

Pressure  R2 

0.991 

0.980 

Figure  24:  Temperature  calibration  curve  for  the  color-camera  method 
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Figure  25:  Temperature  calibration  curve  for  the  two-camera  method 
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Chapter  5:  Experimental  Setup 


Section  5.1:  Experimental  Setup  Introduction 

Three  different  experimental  setups  were  employed  in  this  work,  the  initial 
characterization  of  the  wind  tunnel,  model  surface  pressure-  and  temperature-sensitive 
paint  measurements,  and  off-body  particle  image  velocimetry  measurements. 

The  same  equipment  was  used  to  measure  tunnel  conditions  for  all  three  setups. 
The  wind  tunnel  is  equipped  with  two  350  kPa  absolute  pressure  transducers  which  are 
attached  to  a  total  and  a  static  pressure  port.  Both  have  a  calibration  uncertainty  of  ±  690 
Pa  and  a  frequency  response  greater  than  300  Hz.  The  total  pressure  probe  located  in 
center  of  the  stagnation  tank  is  connected  to  the  pressure  transducer  via  a  0.75-m-long 
6.4mm-diameter  tube.  This  pressure  line  adds  a  phase  lag  to  the  total  pressure  signal 
which  distorts  the  signal  at  frequencies  above  5  Hz.  For  this  reason,  the  total  pressure 
measurement  was  phase  shifted  to  align  with  the  static  pressure  in  order  to  calculate  the 
tunnel’s  Mach  number.  A  Type  K  thennocouple  located  in  the  stagnation  tank  was  used 
to  measure  the  total  temperature.  Isentropic  relations  were  then  used  to  calculate  the  static 
temperature  and  sonic  velocity  in  the  test  section. 
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Section  5.2:  Pressure  and  Temperature  Sensitive  Paint  Experimental 

Setup 


The  airfoil  is  a  NACA  0012  with  a  chord  of  12.7  cm  (5  in)  and  a  15  cm  (6  in) 
span,  as  seen  in  Figure  26.  It  is  equipped  with  10  pressure  taps  of  1.27  mm  inner  diameter 
along  the  upper  surface  and  two  on  the  lower  surface  at  coordinates  listed  in  Table  2. 
The  tap  locations  were  selected  by  the  author  specifically  for  PSP  calibration.  These  taps 
were  connected  to  a  Pressure  Systems,  Inc.  9016  pressure  brick  of  16  differential  pressure 
transducers  with  a  range  of  ±  310  KPa,  an  accuracy  of  ±  210  Pa,  a  measured  frequency 
response  greater  than  50  Hz  and  a  phase  delay  of  7  ms  due  to  acoustic  effects  caused  by 
25.4  cm  of  1.27  mm  inner  diameter  tubing.  These  pressure  taps  were  used  to  perform  an 
in  situ  calibration  of  the  PSP.  Additionally,  two  Te  Technology,  Inc.  MP-2444 
thennistors  were  located  at  x/c  =  0.20,  0.85  and  z^  =  0.933  and  0.933,  respectively,  and 
were  similarly  used  for  an  in  situ  calibration  of  the  TSP.  The  thennistors  were  mounted 
flush  to  the  surface  of  the  airfoil,  had  a  diameter  of  0.89  mm  and  a  range  of  -20  to  100°C. 
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Figure  26:  Diagram  of  the  NACA  0012  airfoil  model  tested 


Table  2:  Pressure  tap  locations  on  the  NACA  0012  airfoil 


Up 

per  Surface 

x/c  = 

0.050 

0.100 

0.150 

0.250 

0.300 

0.350 

0.400 

0.500 

0.700 

0.900 

z/b  = 

0.583 

0.583 

0.583 

0.500 

0.583 

0.667 

0.750 

0.583 

0.583 

0.583 

Lower  Surface 

x/c  = 

0.100 

0.500 

z/b  = 

0.583 

0.583 

A  500-counts-per-revolution  optical  encoder  was  utilized  to  phase  lock  the  PSP 
images  to  the  vane  oscillation.  This  was  accomplished  through  the  use  of  a  National 
Instruments  PCI-6602  counter  board.  The  counter  board  read  the  once-  and  500-per- 
revolution  signals  and  then  sent  a  variable  phase  delay  trigger  pulse  to  both  cameras.  This 
was  done  to  ensure  accurate  phase  separation  for  the  images  while  minimizing  the  run 
duration.  The  total  temperature,  TTL  trigger  pulses,  once  per  revolution,  500  per 
revolution,  thermistors,  total  and  static  pressures  data  were  then  recorded  by  a  National 
Instruments  PCI-625 1  multifunction  data  acquisition  board  to  a  computer  at  a  frequency 
of  4  kHz  for  post  processing. 


48 


A  fast-response,  two-color  porous  polymer  paint  was  selected  for  this  experiment, 
which  was  recently  developed  by  Peng  et  ah. 47  This  paint  is  a  polymer/ceramic  based 
PSP,  which  was  designed  specifically  to  improve  the  diffusivity  of  the  binder  material, 
thus  allowing  faster  response  times.  PtTFPP,  a  typical  platinum  porphyrin  based 
luminophore,  was  used  as  the  pressure  sensor  in  this  experiment.  Perylene  and  a  laser  dye 
were  added  as  the  temperature  sensor.  The  laser  dye  shifts  the  temperature  signal  to  a 
shorter  wavelength  for  improved  spectral  separation  between  the  PSP  and  TSP,  which 
allows  for  the  use  of  optical  filters  to  split  the  data  acquisition  into  red  and  green  signals. 
The  intensity-based  method  utilizes  a  continuous  illumination  source  of  405-nm- 
wavelength  light-emitting  diodes.  The  signals  from  pressure  and  temperature  channels 
are  well  separated  in  spectrum  and  can  be  recorded  simultaneously  either  by  a  color- 
camera  or  a  two-camera  system.  By  co-locating  fast  responding  pressure  and  temperature 
sensors  it  becomes  possible  to  make  high  fidelity  spatial  measurements  of  the  airfoil 
surface  pressure  and  temperature. 

The  two  camera  technique  uses  a  beam  splitter  to  allow  for  both  the  temperature 
and  pressure  cameras  to  have  the  same  optical  view  through  the  limited  optical  access  on 
top  of  the  tunnel.  Each  camera  utilized  an  optical  filter  to  record  the  red  and  green 
illumination  intensities.  The  405  mn  excitation  light  was  provided  by  LED  arrays  through 
the  side  windows  of  the  wind  tunnel.  The  schematic  diagram  of  the  two  camera  setup  can 
be  seen  in  Figure  27.  The  red  channel  utilized  a  590  mn  long-pass  filter  while  the  green 
channel  used  a  5 15  -  600  mn  band-pass  filter  to  separate  the  signals. 
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Figure  27:  Block  diagram  of  the  PSP  data  acquisition  setup,  viewing  down  the 
longitudinal  axis  of  the  wind  tunnel. 


Section  5.3:  Particle  Image  Velocimetry  Experimental  Setup 

Two-dimensional  velocity  data  were  acquired  for  the  NACA  0012  airfoil  using  a 
single  camera  positioned  beside  the  test  section  by  Gompertz.  The  PIV  setup  and  data 
collection  was  performed  by  Kyle  Gompertz,  not  the  author,  and  is  used  as  a  comparison 
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to  the  PSP  and  TSP  data.  The  camera  field  of  view  extends  from  approximately  15%  x/c 
upstream  of  the  airfoil  leading  edge  to  30%  x/c,  to  observe  transient  shock-boundary 
layer  interaction.  The  particle  image  velocimetry  (PIV)  illumination  was  provided  by  a 
200  mJ  /  pulse,  532  nm  wavelength  double-pulsed  Nd:YAG  laser  with  a  maximum  rate 
of  15  Hz.  The  laser  beam  was  expanded  into  a  light  sheet  through  a  cylindrical  negative 
lens  and  entered  the  test  section  from  the  top  window.  The  camera  imaged  the  laser 
sheets  through  the  side  window  of  the  tunnel  test  section  with  a  band-pass  filter  to  only 
record  the  laser  light. 

The  two  consecutive  laser  sheets  were  adjusted  and  measured  to  be  1mm  thick  in 
the  x-y  plane  at  the  airfoil  midspan  location,  z/b  =  0.5.  The  time  separation  between  laser 
sheet  pulses  was  varied  from  0.8  -  1.2  ps  to  limit  the  maximum  particle  displacement  to  7 
pixels  depending  on  the  maximum  time-averaged  velocity  in  the  data  window.  The  flow 
was  seeded  with  Di-Ethyl-Hexyl-Sebacate  (DEHS)  particles  of  submicron  diameter.  As 
described  in  Section  5.2,  an  optical  encoder  on  the  choke  vane  was  used  to  trigger  phase- 
locked  data  acquisition  at  12  discrete  phases  of  the  Mach  number  oscillation  period.  Data 
was  recorded  at/=  9.5  Hz,  a  =  9  and  10°.  Approximately  100  PIV  image  pairs  were 
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acquired  for  each  phase,  which  is  sufficient  to  represent  the  phase-averaged  velocity.  ’ 
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Chapter  6:  Data  Processing  Techniques 


Section  6.1:  Data  Acquisition  Board  Post-Processing  Techniques 

The  information  recorded  from  the  data  acquisition  board  was  saved  to  a 
computer  hard  drive  for  post-processing.  The  first  step  of  the  process  is  to  low  pass  filter 
the  pressure  and  temperature  data.  These  values  along  with  the  angular  encoder  signal 
were  then  averaged  over  the  same  duration  as  the  camera  shutters  were  open.  This 
effectively  created  point  measurements  for  which  to  check  the  a  priori  calibration  of  the 
PSP  and  TSP  and  use  for  an  in  situ  calibration.  In  addition,  these  data  were  also  used  in 
calculating  the  traditional  aerodynamic  condition  data:  Mach  number,  Reynolds  number, 
q ,  P0,  Ps,  and  T0.  The  PSP  and  TSP  images  are  registered,  median  filtered  and  calibrated 
as  discussed  in  Sections  6.2  through  6.4  respectively. 


Section  6.2:  Image  Registration  Techniques 

The  PSP  and  TSP  images  must  be  registered  to  the  same  field  in  order  to  take  the 
appropriate  ratios  and  obtain  reliable  data.  The  wind-on  and  wind-off  images  for  both  the 
green  and  red  channels  must  be  aligned  in  this  way.  Additionally  all  of  the  images  in  the 
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run  must  be  properly  aligned  to  obtain  meaningful  measurements  throughout  the  run. 
This  was  accomplished  by  first  locating  the  fiducial  marker  locations.  Due  to  the  large 
number  of  images,  200  to  1000  per  run,  these  steps  were  automated  and  the  code  can  be 
seen  in  Appendix  B.  The  fiducial  markers,  including  the  pressure  taps  and  airfoil  leading 
and  trailing  edges  were  located  by  taking  the  local  standard  deviation  of  each  wind-on 
and  off  image,  as  seen  in  Figure  28.  A  two-dimensional  cross  correlation  was  performed 
between  the  first  wind-off  standard  deviation  image  and  every  other  image  in  the  run. 
This  forced  all  of  the  images  for  each  run  to  align  with  one  another.  The  images  were 
then  rotated  to  a  right  angle  and  cropped  to  only  show  the  airfoil  surface. 


Figure  28:  Representative  standard  deviation  of  a  wind-off  image  (axes  in  pixels), 
locating  the  fiducial  marker,  pressure  taps  and  airfoil  edges 
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Section  6.3:  Image  Filtering  Techniques  and  Implementation 


Filtering  is  a  typical  data  processing  technique  in  pressure  sensitive  paint; 
however,  there  are  tradeoffs  between  different  filter  types  which  will  be  described  in  this 
section.  By  applying  a  top  hat  filter  over  the  image  surface  one  is  essentially  low-pass 
filtering  the  data  with  the  goal  of  removing  the  erroneous  white  noise  while  maintaining 
the  data  quality.  To  illustrate  this  process  a  nominal  PSP  image  was  selected  from  a 
steady  run  at  M  =  0.46,  Re  =  3.65><106  and  a  =  9°,  as  seen  in  Figure  29.  This  image  was 
chosen  due  to  the  relatively  weak  standing  shock  on  the  leading  edge  of  the  airfoil. 
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Figure  29:  Unfiltered  Cp  map  of  the  upper  surface  NACA  0012  at  M  =  0.46,  Re  = 
3.65x106  and  a  =  9°,  flow  from  left  to  right 


Three  standard  digital  filters  were  applied  in  the  post-processing  of  the  PSP 
images  to  evaluate  the  performance  of  each  filter,  mean,  Wiener,  and  median.  The  mean 
filter  as  seen  in  Figure  30  used  a  disk  top  hat  filter  with  a  radius  of  8  pixels.  The  Wiener 
and  median  filters  seen  as  Figure  31  and  Figure  32  respectively  used  a  rectangular  filter 
pattern  of  30  pixels  in  the  spanwise  direction  and  10  pixels  in  the  chordwise  direction. 
The  filter  was  shorter  in  the  chordwise  direction  due  to  the  higher  changes  in  pressure  as 
compared  to  that  of  the  spanwise  direction. 
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Figure  30:  Mean  top  hat  filtered,  with  a  radius  of  8  pixels,  Cp  map  of  the  upper  surface 
NACA  0012  at  M  =  0.46,  Re  =  3.65xl06  and  a  =  9° 
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Figure  3 1 :  Wiener  top  hat  filtered,  30  pixels  in  the  span  wise  and  10  pixels  in  the  chord 
wise  direction,  Cp  map  of  the  upper  surface  NACA  0012  at  M  =  0.46,  Re  =  3.65x  106  and 
a  =  9° 
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Figure  32:  Median  top  hat  filtered,  30  pixels  in  the  span  wise  and  10  pixels  in  the  chord 
wise  direction,  Cp  map  of  the  upper  surface  NACA  0012  at  M  =  0.46,  Re  =  3.65x  106  and 
a  =  9° 


These  raw  and  filtered  images  were  then  compared  directly  for  one  Cp  slice  at  the 
spanwise  location  of  350  pixels  or  z/b  =  0.608,  just  below  the  majority  of  the  pressure 
taps  at  z/b  =  0.583.  Figure  33  shows  the  full  chord  for  all  of  the  filters  while  Figure  34  is 
zoomed  in  on  the  first  25%  of  the  NACA  0012.  It  becomes  apparent  that  all  three  filters 
remove  a  majority  of  the  high  frequency  noise.  However,  the  median  filter  preserves  the 
shock  gradient  at  the  leading  edge  and  reduces  the  white  noise  most  effectively. 
Furthermore,  median  filters  are  standard  practice  for  non-parametric  data  analysis  where 
a  normal  distribution  cannot  be  assumed.49  Thus  with  varying  sample  sizes,  an  unknown 
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normal  distribution,  and  the  ability  to  preserve  abrupt  features  such  as  shocks  while 
eliminating  high  frequency  noise,  the  median  filter  was  chosen  as  the  image  processing 
filter  for  this  experimental  data  set.  With  the  setup  used  there  were  5.14  pixels  per 
percent  x/c. 


-4 
-3.5 
-3 
-2.5 
-2 

0-15 
-1 
-0.5 
0 
0.5 
1 

0  100  200  300  400  500 

Pixels 


No  Filter 
Mean  Filter 
Wiener  Filter 


Median  Filter 


Figure  33:  Cp  profile  plot  of  the  unfiltered  and  three  filtering  methods  of  a  NACA  0012  at 
M=0.46,  Re=3.65><106  and  a=9°  at  the  spanwise  location  of  350  pixels 
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Figure  34:  First  25%  x/c  Cp  profile  plot  of  the  unfiltered  and  three  filtering  methods  of  a 
NACA  0012  at  M=0.46,  Re=3.65x106  and  a=9°  at  the  spanwise  location  of  350  pixels 


Section  6.4:  Pressure  and  Temperature  Calibration 

A  priori  and  in  situ  calibrations  are  the  two  major  methods  of  calibration  for 
pressure  sensitive  paint.  The  a  priori  method  utilizes  a  pressure  and  temperature 
calibration  chamber  to  record  the  emission  intensity  of  the  PSP  and  TSP  at  known,  steady 
pressures  and  temperatures.  These  intensities  are  then  divided  by  the  reference  condition 
to  calculate  the  intensity  ratio  for  each  point  in  pressure  and  temperature.  This  method 
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was  attempted  for  this  experimental  setup,  although  changes  in  the  reference  condition 
carried  large  impacts  into  the  results.  The  reference  condition  in  the  a  priori  method  is 
assumed  to  be  at  a  uniform  pressure  and  temperature,  which  is  not  valid  after  even  one 
run.  As  cold,  expanded  air  flows  through  the  tunnel,  the  airfoil  and  2.54-cm-thick 
aluminum  side  walls  cool  which  is  slowly  dissipated  back  to  the  atmosphere.  This  causes 
the  airfoil  to  have  a  temperature  gradient  during  the  wind-off  image  as  the  thin  trailing 
edge  cools  and  warms  more  rapidly  than  the  thicker  leading  edge.  Additionally  the  sides 
of  the  airfoil  model  are  different  as  seen  in  Figure  26,  which  causes  a  spanwise  gradient 
in  temperature.  This  heat  transfer  between  runs  negates  the  assumption  of  uniform 
temperature  over  the  reference  image. 

For  a  bi-luminophore  PSP,  the  in  situ  method  is  implemented  by  defining  the 
calibration  of  PSP  intensity  ratio  based  on  the  surface  pressure  taps.  However  a  typical  in 
situ  calibration  was  not  possible  in  this  test  due  to  the  attenuation  and  phase  shift  in  the 
pressure  tap  lines  during  the  unsteady  Mach  oscillation  runs.  For  this  reason  the  PSP  data 
could  not  be  calibrated  for  the  high  frequency  Mach  oscillation  runs  and  was  only  used  to 
investigate  the  shock  location. 

A  surface  fit  calibration  method  was  devised  to  compensate  for  the  transient 

temperature  field.  This  in  situ  method  utilizes  pressure  tap  data,  of  the  steady  runs,  to 

create  a  surface  fit  calibration  for  pressure  with  respect  to  the  green  and  red  channel 

intensity  ratios.  This  surface  fit  calibration  method  allows  for  one  calibration  for  each  run 

that  accounts  for  both  variations  in  pressure  and  temperature.  Whereas,  with  single- 

luminophore  pressure-sensitive  paints  the  calibration  changes  with  variations  in 

temperature,  requiring  a  new  fit  for  each  image.  An  example  surface  fit  calibration  can  be 
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seen  in  Figure  35.  While  Figure  36  shows  how  the  calibration  is  implemented  to  remove 
the  temperature  sensitivity  from  the  PSP  measurement.  This  calibration  is  implemented  to 
find  the  coefficient  of  pressure  for  a  representative  steady  run  as  seen  in  Figure  37. 
Historic  data  is  also  presented  for  comparison  and  is  discussed  fully  in  Section  7.1. 
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Figure  35:  PSP  in  situ  surface  fit  calibration  for  pressure  vs.  green  and  red  intensity  ratios 
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Figure  36:  Schematic  diagram  of  bi-luminophore  pressure-  and  temperature-sensitive 
paint  in  situ  pressure  calibration 
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Figure  37:  Three  Cp  profiles  from  calibrated  PSP  data  of  the  NACA  0012  airfoil  at  M  = 
0.6,  a  =  9°  and  Re  =  3.8  million  shown  with  a  Cp  profile  of  a  NACA  0012  M  =  0.6,  a  = 
9°  and  Re  =  3  million 


Section  6.5:  Experimental  Uncertainty 

Several  factors  contribute  to  the  PSP  experimental  uncertainty  such  as  the 
pressure  tap  accuracy,  model  movement  and  defonnations.  PSP  is  also  susceptible  to 
variations  in  the  illumination  field  which  can  occur  from  model  movement  through  a 
non-unifonn  excitation  light.  This  is  minimized  by  rigidly  affixing  the  cameras,  beam 
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splitter,  airfoil  and  LED  arrays  to  the  tunnel.  Another  source  of  illumination  error  is 
caused  by  variations  in  broad  spectral  noise  from  changes  in  ambient  illumination 
intensities;  these  unwanted  background  illuminations  are  minimized  by  limiting  the  light 
entering  through  the  tunnel  windows  by  turning  off  the  room  lights  or  covering  the 
window  with  a  blackout  curtain.  However,  due  to  the  tunnel  exhausting  to  the  outdoors, 
some  daylight  enters  the  tunnel.  This  background  illumination  field,  typically  called  the 
dark  image,  is  recorded  and  subtracted  from  the  wind-on  and  wind-off  images.  Only  one 
set  of  background  images  is  taken  for  each  run  and  it  is  assumed  that  the  ambient  light 
does  not  change  from  the  wind-on  to  wind-off  images.  Image  registration  between  the 
green  wind-on,  green  wind-off,  red  wind-on  and  red  wind-off  images  is  also  a  source  of 
error  due  to  the  non-uniform  pressure,  temperature,  luminophore  thickness,  and 
illumination.  By  incorrectly  aligning  all  of  the  images  these  error  effects  can  be 
exacerbated. 


Section  6.6:  Shock  Location  Measurement  Techniques 


PSP,  TSP  and  PIV  images  have  been  acquired  for  different  azimuthal  locations, 
angles  of  attack  (a),  and  Mach  oscillation  frequencies.  Each  image  was  interrogated  to 
find  the  location  of  the  shock.  With  the  high  spatial  resolution  of  the  pressure  sensitive 


paint  the  shock  location  is  found  by  taking  the  derivative  of  Cp  with  respect  to  x/c, 


dCp 
dx/c ’ 


fiT 

and  likewise  the  derivative  of  temperature  with  respect  to  x/c,  ~~r-  The  shock  location 
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for  PSP,  PIV  and  TSP  is  thus  found  by  locating  the  derivative  extremum  within  the 
bounded  neighborhood  where  a  shock  is  expected,  3  <  x/c  <  30%,  as  expressed  in 

17 

Equation  9  and  visualized  in  Figure  38.  The  reported  shock  position  is  the  location 
where  the  shock  strength  was  found  to  be  at  a  minimum,  though  multiple  pixels  were 
found  to  be  within  the  shock.  In  each  case,  a  peak  in  the  derivative  is  identified  as  the 
location  of  the  shock.  The  automated  data  reduction  techniques  necessary  to  obtain 
meaningful  shock  location  results  for  large  data  sets  are  presented  in  Appendices  A 
through  E. 


Equation  9:  Shock  location  equation. 
Shock  Location  =  min  { 

d  x  c 
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Figure  38:  Example  upper  surface  Cp  profile  and  derivative  of  the  Cp  with  respect  to  x/c 

Figure  39  shows  a  representative  suction-surface  Cp  distribution  acquired  using 
the  single-shot,  two-camera,  temperature  corrected,  intensity-based  PSP  method;  the 
shock  location  can  be  seen  as  the  abrupt  jump  in  Cp  between  -4  and  -3  at  x/c  =  13%.  The 
shock  location  of  the  PSP  and  TSP  images  were  analyzed  along  the  centerline  location,  in 
order  to  match  the  PIV  laser  plane.  Figure  39  shows  a  representative  surface  contour  Cp 
map  with  the  interrogation  location  seen  as  a  black  line.  Figure  40  shows  one  phase  of 
ensemble-averaged  phase-locked  PIV  images.  Note  the  PSP  and  PIV  representative 
images  are  at  two  different  test  conditions.  The  black  line  marks  the  line  of  interrogation 
for  the  shock  location;  this  measurement  is  made  approximately  4  mm  above  the  airfoil 
due  to  reflected  laser  light  off  the  airfoil,  giving  erroneous  vector  fields  close  to  the 

T9 

surface. 
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Figure  39:  A  representative  suction  surface  Cp  distributions  of  a  NACA  0012  at  a  =  9°,  M 
=  0.55  and  Re  =  3.6  million,  flow  is  from  left  to  right.  The  black  line  is  where  the  shock 
location  was  searched  for 
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Figure  40:  Contour  of  ensemble-averaged  Mach  number  for  one  phase  of  a  Mach 
oscillation  period  with  M  =  0.6,  a  =  9°  and /=  9.5  Hz.  The  black  line  is  where  the  shock 
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location  was  searched  for 
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Chapter  7:  Results 


Section  7.1:  Static  Results 

The  NACA  0012  airfoil  was  tested  at  9,  10  and  1 1  degrees  angle  of  attack  and  at 
Mach  numbers  of  0.46,  0.55  and  0.61.  These  Mach  numbers  were  selected  to  bound  the 
unsteady  data  in  order  to  make  direct  comparisons,  while  also  at  the  same  angles  of 
attack.  Additionally  the  M  =  0.61  data  at  a  =  10°  is  in  the  buffeting  regime  and  at  an 
angle  of  attack  of  =  11°  the  airfoil  is  stalled.  The  corresponding  Reynolds  numbers 
varied  between  3.4  and  4.4  million  for  the  0.127  m  chord.  The  measurements  from  the  bi- 
luminophore  pressure-  and  temperature-sensitive  paint  were  in  situ  calibrated  to  pressure 
and  then  converted  into  a  coefficient  of  pressure  using  Equation  10.  The  Cp  values  were 
then  converted  to  a  normal  ( C„)  and  axial  force  coefficients  (Ca)  by  integrating  along  the 
airfoil  chord  as  seen  in  Equation  1 1  and  Equation  12.  Further,  the  data  was  converted  into 
sectional  coefficient  of  lift  (C/),  drag  (CJ)  and  moment  (Cm)  by  means  of  Equation  13 
through  Equation  15.50  These  calculations  were  perfonned  only  using  the  pressure  forces 
and  neglecting  the  friction  forces  on  the  airfoil.  The  results  from  these  calculations  are 
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shown  in  Appendix  F,  Table  4.  ’ 
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Equation  10:  Coefficient  of  pressure 


C 


P 


f  \ 


yps  ) 


1/ 
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yM2 


Equation  1 1 :  Normal  force  coefficient  equation 
1 

c„  =  |  {cpj  -  Cp  u  )d  x/c 
0 


Equation  12:  Axial  force  coefficient  equation 


dyu/c 

dx/c 


C 


pi 


dyjc 
d  x/c 


\ 

d  x/c 


Equation  13:  Sectional  coefficient  of  lift  equation 
C,  =  Cn  cos ci-Ca  sina 


Equation  14:  Sectional  coefficient  of  drag  equation 
Cd  =  Cn  sin  a  +  Ca  cos  a 
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Equation  15:  Sectional  coefficient  of  moment  equation  about  the  quarter  chord 

0.25  1 

c,  =  I  (c„  -  C„  Xo.25  -  x/c>  x/c  -  J  (c„  -  C„  )(i/c -0.25)rf  x/c 

0  0.25 


Three  representative  Cp  profiles  are  shown,  Figure  41  through  Figure  48,  for  each 
static  test  condition  representing  the  data  acquired  at  10,  50,  and  90%  of  the  run  duration. 
This  effectively  shows  the  Cp  profile’s  insensitivity  towards  temperature,  as  the  airfoil  is 
cooled  throughout  the  run.  Pressure  taps  for  all  three  Cp  profiles  are  shown  in  the  same 
color  as  the  PSP  data.  The  Cp  data  after  the  shock  on  the  upper  or  suction  surface  should 
be  trending  in  a  straight  line  towards  Cp  =  0  as  seen  with  the  historical  NASA  data  in 
Figure  37. 
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Figure  41 :  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  9°,  M  = 
0.46  and  Re  =  3.5  million  with  pressure  taps  shown 


Figure  42:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  10°,  M 
=  0.45  and  Re  =  3.4  million  with  pressure  taps  shown 
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Figure  43:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  1 1°,  M 
=  0.45  and  Re  =  3.4  million  with  pressure  taps  shown 


Figure  44:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  9°,  M  = 
0.55  and  Re  =  3.7  million  with  pressure  taps  shown 
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Figure  45:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  10°,  M 
=  0.55  and  Re  =  3.7  million  with  pressure  taps  shown 


Figure  46:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  1 1°,  M 
=  0.55  and  Re  =  3.8  million  with  pressure  taps  shown 
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Figure  47:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  9°,  M  = 
0.61  and  Re  =  4.1  million  with  pressure  taps  shown 


Figure  48:  Three  representative  PSP  Cp  profiles  of  the  NACA  0012  airfoil  at  a  =  10°,  M 
=  0.61  and  Re  =  4.3  million  with  pressure  taps  shown 
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The  Ci,  Cm,  C„  and  Ca  were  calculated  for  all  50  to  100  image  pairs  depending  on 
the  run.  The  standard  deviation  was  then  calculated  for  each  set  of  data,  in  order  to 
evaluate  the  precision  error  for  each  measurement  technique.  Similarly  a  direct 
comparison  of  the  M  =  0.6  PSP  and  NASA  data  is  performed  to  evaluate  the  bias  error  in 
the  measurement.  At  M  =  0.6  and  a  =  9°,  the  PSP  results  give  a  C/  =  0.882  ±  0.022  and 
Cm  =  0.015  ±  0.012  while  the  NASA  data  gives  C/  =  0.852  and  Cm  =  0.0232.  This  results 
in  a  bias  error  of  3.5%  and  89%  for  the  Ci  and  Cm  measurements  respectively.  The 
Coefficient  of  lift  and  moment  vs.  Mach  number  can  be  seen  in  Figure  50  and  Figure  51 
for  all  static  test  conditions,  the  raw  data  can  be  found  in  Appendix  F  as  Table  4.  It  is 
expected  that  the  coefficient  of  moment  would  have  much  higher  errors  than  the 
coefficient  of  lift  PSP  measurements  due  to  the  error  in  trailing  edge  from  lingering 
temperature  effects  and  low  SNR  from  ambient  light.  These  effects  can  be  seen  in  the  Cp 
plots  as  waviness  at  x/c  >  70%.  This  waviness  is  due  to  the  low  illumination  field  at  the 
trailing  edge  of  the  airfoil,  as  shown  in  Figure  49  where  the  ratio  of  background  light  to 
wind-off  illumination  reaches  14%  at  the  trailing  edge  of  the  airfoil.  Additionally  the  first 
few  percent  of  the  leading  edge  is  difficult  to  measure  measure/image  for  the  negative 
angles  of  attack.  The  NASA  Technical  Memorandum  data  utilized  a  wake  probe  to 
calculate  drag  while  the  PSP  data  only  accounts  for  pressure  drag  and  not  skin  friction. 
Therefore  it  can  be  expected  that  the  NASA  and  PSP  values  of  O  would  differ 
substantially  and  as  such  are  not  presented. 
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Pixels 


Pixels 

Figure  49:  Dark  image  divided  by  wind-off,  showing  the  large  susceptibility  to  varying 
background  illumination  at  the  trailing  edge 
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Figure  50:  Coefficient  of  lift  vs.  Mach  number  for  a  =  9,  10  and  1 1°,  for  both  PSP  and 
historical  data 
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Figure  5 1 :  Coefficient  of  moment  vs.  Mach  number  for  a  =  9,  10  and  1 1°°,  for  both  PSP 
and  historical  data 


Section  7.2:  Buffeting  Results 

The  NACA  0012  airfoil  is  within  the  buffeting  regime  at  the  static  tested 
condition  of  a  =  10°  and  M  =  0.61,  as  seen  earlier  in  Figure  48.  This  critical  test 
condition  demonstrates  the  utility  of  the  bi-luminophore  PSP  to  detect  and  track  this 
highly  unsteady,  three-dimensional  phenomenon.  Figure  52  shows  9  representative  upper 
surface  Cp  maps  for  0%  <  x/c  <  60%  and  20%  <  z/b  <  80%.  These  successive  images,  left 
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to  right  and  then  top  to  bottom,  show  the  highly  3D  nature  of  the  flow.  These  bi- 
luminophore  images  had  an  exposure  time  of  5.5  ms  for  both  the  green  and  red  channels 
and  were  acquired  at  a  rate  of  28  Hz.  While  the  PSP’s  frequency  response  has  been 
characterized  to  be  greater  than  5  kHz,  the  technique  is  limited  by  the  duration  needed  to 
adequately  record  the  emission  intensity  signals.  The  single-shot  unsteady  PSP  method  is 
pivotal  to  measuring  this  highly  dynamic  phenomenon.  Any  non-phase-locked  PSP 
techniques  would  have  averaged  out  this  fluid  dynamic  phenomenon,  if  the  single-shot 
technique  was  not  utilized. 
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Figure  52:  9  upper  surface  Cp  maps  of  the  NACA  0012  airfoil  at  a  =  10°,  M  =  0.61  and 
Re  =  4.2  million  for  0%  <  x/c  <  60%  and  20%  <  z/b  <  80%  during  a  steady  run  in  the 
buffeting  regime 
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Section  7.3:  Dynamic  Results 


In  addition  to  the  static  testing,  the  NACA  0012  airfoil  was  tested  at  a  variety  of 
low-reduced-frequency  Mach  oscillation  rates  and  angles  of  attack,  as  seen  in  Table  3. 
The  mean  Mach  number  was  utilized  to  find  the  reduced  frequency  for  the  unsteady  data. 
The  objective  of  this  research  was  to  evaluate  the  shock  movement  in  conjunction  with 
the  coefficient  of  lift,  moment  and  drag  hysteresis  curves  for  each  angle  of  attack  and 
frequency.  It  was  not  possible  to  obtain  accurate  pressure  distributions  on  the  airfoil  due 
to  the  lack  of  unsteady  pressure  transducers  for  an  adequate  in  situ  calibration,  in 
conjunction  with  the  poor  quality  of  the  a  priori  calibration.  Therefore,  the  unsteady  data 
is  limited  to  unsteady  shock  location  detennination  between  the  PSP,  TSP  and  PIV  in 
addition  to  qualitative  analysis  of  the  fluid  dynamic  phenomenon. 


Table  3:  Unsteady  test  matrix 


Oscillation  Frequency  (Hz) 

Reduced  Frequency 

Angle  of  Attack  (°) 

21 

0.050 

9,  10  and  11 

15.25 

0.037 

9,  10  and  11 

9.5 

0.023 

9,  10  and  11 

2.1 

0.005 

9,  10  and  11 

The  Mach  oscillation  waveforms  were  correlated  to  azimuthal  angle  by 
perfonning  a  cross  correlation  of  the  data  to  align  with  a  sinusoidal  wavefonn  which  has 
a  maximum  Mach  number  at  W  =  90°.  This  was  done  to  enable  a  direct  comparison 
between  oscillation  frequencies,  due  to  the  varying  phase  shift  with  frequency,  as 
discussed  in  Section  3.2.  The  azimuthally  aligned  data  can  be  seen  in  Figure  53  through 
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Figure  55  corresponding  to  a  =  9,  10  and  11°  respectively,  with  each  data  point 
corresponding  to  the  Mach  number  at  each  PSP  image.  The  maximum  frame  rate  for 
both  cameras  was  28  Hz  for  these  tests,  thus  Figure  53  through  Figure  61  are  a 
compilation  of  multiple  oscillations  periods.  Due  to  the  addition  of  an  airfoil  in  the  tunnel 
the  Mach  oscillation  wavefonn  was  shifted  from  the  empty  test  section  waveforms  (see 
Section  3.2).  This  shift  is  due  to  the  change  in  the  cross  sectional  area  in  the  test  section, 
At.  This  effect  can  be  seen  in  Figure  53  through  Figure  55  as  a  bulk  shift  in  the  Mach 
oscillation  waveforms  due  to  the  change  in  angle  of  attack  from  a  =  9  to  11°.  The  bulk 
shift  was  down  which  is  counter-intuitive  to  the  Mach  -  area  equation,  Equation  1 .  This  is 
possibly  due  to  the  static  pressure  transducer  being  too  close  to  the  airfoil.  Additionally 
with  the  pressure  transducer  being  located  upstream  of  the  airfoil,  the  cross  sectional  area 
is  slightly  increased  thereby  emulated  a  lower  Mach  number. 
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Mach  number  £  Mach  number 


53:  Mach  number  vs.  Phase  for  the  4  frequencies  at  a  =  9° 
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Figure  54:  Mach  number  vs.  Phase  for  the  4  frequencies  at  a  =  10° 
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Figure  55:  Mach  number  vs.  Phase  for  the  4  frequencies  at  a  =  11° 


The  Mach  oscillation’s  effect  on  the  shock  location  was  found  to  have  a 
significant  phase  delay  associated  with  the  frequency.  This  effect  can  be  seen  as  a  growth 
in  hysteresis  in  Figure  56  with  an  increase  in  frequency.  This  phenomenon  was  also 
measured  with  unsteady  Schlieren  images  by  Fernie  and  Babinsky  at  similar  low  reduced 
frequencies  on  the  NACA  0012.  ’  For  the  high  frequency  Mach  oscillations  at  a  =  9°, 
the  retreating  shock  is  found  to  lag  the  advancing  shock,  Figure  56.  The  hysteresis  effect 
is  due  to  the  instantaneous  Mach  number  varying  and  the  fluid  dynamic  phase  lag  of  the 
shock  location.  The  start  of  the  schematic  directional  loops  in  this  figure  correspond  to  ¥ 
=  0°  and  represents  the  advancing  blade  through  ¥=  180°  while  ¥=  180°  through  360° 
correspond  to  the  retreating  blade.  The  a  =  10°  data,  Figure  57,  is  interesting  when 
compared  to  the  steady  data,  for  which  the  M  =  0.61  data  was  found  to  be  buffeting  while 
the  M  =  0.55  was  steady  attached  flow.  It  can  be  seen  that  at  all  frequencies  at  and  above 
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M  =  0.55  the  shock  is  buffeting,  however  as  the  frequency  increases  the  steady-to- 
buffeting  transition  occurs  at  lower  Mach  number.  It  is  hypothesized  that  the  forced 
freestream  flow  causes  the  turbulent  boundary  layer  to  breakdown  into  a  buffeting  flow 
more  rapidly.  Figure  58  shows  a  continuation  of  this  aperiodic  trend  with  even  the  low 
Mach  numbers  exhibiting  cycle-to-cycle  variations  in  the  shock  location.  Figure  59 
through  Figure  61  illustrate  this  phenomenon  as  shock  location  vs.  azimuth  angle,  where 
it  is  seen  that  the  shock  location  is  very  steady  at  the  retreating  blade  angles  around  ¥  = 
270°,  where  the  Mach  number  is  lower. 
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Figure  56:  Mach  number  vs.  Shock  Location  for  the  4  frequencies  at  a  =  9°,  with  hand 
drawn  schematic  directional  loops  of  the  oscillations  starting  at  W=  0° 
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Figure  57:  Mach  number  vs.  Shock  Location  for  the  4  frequencies  at  a  =  10° 


87 


o 

2.1  Hz 

X 

9.5  Hz 

V 

15.25  Hz 

□ 

21  Hz 

Figure  58:  Mach  number  vs.  Shock  Location  for  the  4  frequencies  at  a  =  11° 


Figure  59:  Shock  Location  vs.  Phase  for  the  4  frequencies  at  a  =  9° 
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gure  60:  Shock  Location  vs.  Phase  for  the  4  frequencies  at  a  =  10° 
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Figure  61 :  Shock  Location  vs.  Phase  for  the  4  frequencies  at  a  =  11° 
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The  buffeting  regime  was  further  investigated  for  the  oscillatory  flow  tests  with 
PSP  and  PIV.  Four  phase-locked,  XP=  93°  and /=  9.5  Hz,  PIV  and  PSP  images  are  shown 
at  a  =  9°  in  Figure  62,  which  is  below  the  buffeting  regime  and  in  the  buffeting  regime  at 
a  =  10°  as  seen  in  Figure  63.  The  PSP  and  PIV  images  are  4  representative  images  of  the 
same  tunnel  conditions  and  azimuth  position  taken  during  different  runs  to  show  the 
periodic  or  aperiodic  nature  of  the  flow.  It  is  shown  that  the  a  =  9°  results  are  highly 
repeatable  with  both  the  on-  and  off-body  measurement  techniques.  The  higher,  a  =  10° 
images  show  an  aperiodic  phenomenon  similar  to  that  seen  in  the  PSP  buffeting  results. 
The  2D  PIV  results  exhibit  variations  in  the  shock  strength  and  location  while  the  PSP 
results  show  this  phenomenon  to  be  three  dimensional  in  nature,  similar  to  what  was 
encountered  in  the  steady-flow  buffeting  results. 
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Figure  62:  Four  representative  phase  locked  images,  W=  93°,  of  PIV  and  PSP  results  at 
9.5  Hz,  a  =  9°  and  M  =  0.59  showing  a  steady  shock  location 
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Figure  63:  Four  representative  phase  locked  images,  ! F=  90°,  of  PIV  and  PSP  results  at 
9.5  Hz,  a  =  10°  and  M  =  0.59  showing  an  unsteady  shock  location 
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Due  to  the  length  and  number  of  tunnel  runs  needed  to  obtain  100  images  the  PIV 
data  has  only  been  recorded  for  one  full  set  of  Mach  oscillation  phases  at  one  frequency. 
This  set  of  ensemble-averaged  PIV  images  for  12  phases  of  a  Mach  oscillation  was  taken 
at  a  =  9°  and /=  9.5  Hz  and  can  be  seen  in  Figure  64.  This  shows  the  off-body  forces  of 
the  Mach  oscillation  progressing  from  a  strong  to  weak  shock. 
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Figure  64:  Contours  of  ensemble-averaged  Mach  number  for  12  phases  of  a  Mach 
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oscillation  period  with  a  =  9°  and /=  9.5  Hz 
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Figure  65:  Mach  number  vs.  Shock  location  for  PSP,  TSP  and  PIV  at /=  9.5  Hz  and  a= 9° 


The  spatial  resolution  of  the  PSP  and  TSP  measurements  of  shock  location  was 
0.125%  x/c  while  PIV  had  a  slightly  lower  resolution  of  0.33%  x/c.  The  TSP 
measurement  was  taken  from  the  reference  channel  emission  (Perylene  and  laser  dye 
mixture)  for  the  PSP  temperature  compensation.  A  direct  comparison  of  the  shock 
detection  results  are  shown  in  Figure  65  for  the  one  Mach  oscillation  frequency  tested 
with  all  three  measurement  techniques  (PIV,  PSP,  and  TSP).  It  can  be  seen  that  the  PSP 
and  PIV  both  show  the  same  trend  while  TSP  is  very  erroneous  in  tracking  the 
instantaneous  shock  location;  the  TSP  measures  the  solid  aluminum  airfoil  surface 
temperature  and  not  the  shock  location  which  are  tenuously  linked  due  to  the  effects  of 
heat  conduction.  The  PIV  measurement  is  seen  to  be  shifted  to  a  higher  chordwise 
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location  at  higher  Mach  numbers  due  to  the  measurement  taken  at  4  mm  above  the  airfoil 
surface  along  with  the  shock  curving  back  to  an  aft  position  as  it  increases  in  height. 
Hysteresis  effects  are  particularly  evident  where  shock  strength  is  shown  to  be  dependent 
upon  whether  the  freestream  is  accelerating  or  decelerating.  This  hysteresis  is  observed 
at  a  low  value  of  reduced  frequency  (0.037),  which  is  normally  considered  to  be  quasi- 

■J  1 

steady.  This  is  possible  due  to  the  ability  to  measure  the  shock  location  with  greater 
accuracy,  thereby  detecting  smaller  unsteady  effects. 


95 


Chapter  8:  Conclusions 


Section  8.1:  Summary  of  Experimental  Results 

A  bi-luminophore  unsteady-pressure-  and  temperature-  sensitive  paint  was  used  in 
conjunction  with  a  new  two-camera  technique  to  study  the  steady  and  unsteady 
characteristics  of  shock  motion  on  the  NACA  0012  airfoil.  In  an  initial  study  a  single¬ 
shot,  lifetime -based,  unsteady  pressure-sensitive  paint  method  was  utilized  to  study  the 
shock  motion  of  a  NACA  0021  airfoil.  This  method  effectively  resolved  the  unsteady 
standing  shock  wave  on  the  airfoil.  However,  due  to  temperature  non-uniformity  in  the 
wind-off  images  and  temporal  gradients  in  the  wind-on  images,  the  trailing  edge  Cp 
profile  could  not  accurately  be  resolved.  For  this  reason  the  single-shot  intensity-based 
technique  with  a  bi-luminophore  fast  pressure-paint  was  investigated  for  use  in  this 
tunnel.  The  high  spatial  resolution  of  PSP,  TSP  and  PIV  measurements  allowed  for  high 
accuracy  of  shock  location  capture.  The  single-shot  PSP  technique  allows  high  accuracy 
in  the  chordwise  and  spanwise  direction  allowing  for  the  capture  of  aperiodic  3D  fluid 
dynamic  phenomenon,  such  as  buffeting.  If  a  phase  averaged  technique  were 
implemented  these  aperiodic  fluid  dynamic  features  of  the  flow  would  have  been  lost  to 
cyclic  averaging.  It  was  found  that  these  perturbations  in  shock  location  are  not  2D  for  a 
2D  airfoil  section.  Bi-luminophore  pressure-sensitive  paint  has  demonstrated  the 
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capability  to  produce  high-resolution  data  with  the  new  two-camera  single-shot  intensity- 
based  technique.  Standard  aerodynamic  forces  (C„,  Ca_  Ci  and  Cm)  were  computed  for 
only  the  steady  runs,  due  to  the  lack  of  adequate  calibration  techniques  for  the  unsteady 
results.  Techniques  and  practices  to  alleviate  these  errors  are  discussed  in  Section  8.2. 
The  steady  forces  were  compared  to  similar  NACA  0012  historical  data.  It  was  found 
that  the  PSP  was  able  to  measure  the  sectional  coefficient  of  lift  with  high  accuracy  and 
precision.  The  errors  of  the  PSP  measurements  were  found  to  be  3.5%  and  89%  for  the  C/ 
and  Cm  measurements,  respectively.  Much  higher  errors  were  found  with  the  coefficient 
of  moment  measurement  due  to  the  errors  in  the  trailing  and  leading  edges.  The 
erroneous  pressure  forces  at  the  leading  and  trailing  edges  were  exacerbated  by  their 
large  moment  arm  relative  to  the  quarter  chord  location.  Additionally,  the  error  seems 
large  due  to  the  small  values  of  Cm  for  a  symmetric  airfoil. 

The  airfoil  was  tested  in  the  buffeting  regime  at  one  steady  condition.  The  shutter 
of  the  camera  was  open  for  approximately  one  third  of  the  buffeting  period  as  compared 
to  historical  buffeting  studies.  Additionally  the  buffet  frequency  was  higher  than  the 
Nyquist  frequency  of  the  pressure  taps.  In  order  to  fully  study  the  three  dimensionality  of 
this  phenomenon  an  optimized  in-depth  study  is  needed.  This  would  consist  of  shorter 
duration  camera  shutters  through  the  use  of  more  excitation  illumination  and  a  high 
frequency  transducers  surface  mounted  on  the  airfoil  for  accurate  in  situ  pressure 
calibration.  However,  the  bi-luminophore  PSP  single-shot  technique  is  uniquely  suited  to 
make  these  3D  aperiodic  measurements. 


97 


Section  8.2:  Recommendations  for  Obtaining  Accurate  Unsteady  Cp 

Data  Using  PSP 

Several  experimental  setup  techniques  were  employed  to  minimize  error. 
Theoretically  is  it  possible  to  perform  and  implement  an  a  priori  calibration  to  obtain  the 
surface  pressure  data;  however,  in  practice  this  calibration  is  done  on  a  coupon  and  not 
the  tested  airfoil.  This  leads  to  variations  in  the  illumination  field,  luminophore 
concentration  and  camera  settings  between  the  a  priori  calibration  coupon  and  the  tested 
airfoil.  These  effects  are  minimized  by  taking  the  ratio  of  the  wind-off  and  wind-on 
images,  but  they  are  not  eliminated  entirely  due  to  registration  errors  or  model 
deformation  and  introduce  errors  into  the  measurement  technique.  Therefore,  in  practice 
in  situ  calibrations  are  almost  always  used,  which  require  accurate  pressure 
measurements  at  all  testing  conditions.  It  is  due  to  the  lack  of  unsteady  surface  mounted 
pressure  transducers  that  unsteady  coefficient  of  lift  and  moment  plots  cannot  be 
presented  in  this  work.  For  bi-luminophore  paint  the  best  calibration  will  arise  from  the 
thennistors  and  pressure  taps  being  spread  in  the  temperature  and  pressure  regime  tested, 
thus  minimizing  any  extrapolation  errors  in  the  calibration  implementation.  It  was  found 
that  the  intensity  of  the  ambient  light  changed  between  the  wind-on  and  wind-off  images 
due  to  variations  in  the  light  entering  the  exhaust  of  the  wind  tunnel.  The  changing  dark 
or  ambient  illumination  field  could  be  accounted  for  by  using  a  double-shutter  CCD 
camera  in  conjunction  with  the  intensity-based  PSP  measurement,  using  Gate  1  for  the 
intensity-based  PSP  measurement  and  Gate  2  as  a  dark  reference  image.  The  dark  image 
would  thus  be  separated  from  the  wind-on  or  wind-off  image  on  the  order  of  a 

millisecond  instead  of  tens  of  seconds  with  the  traditional  technique. 
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The  two-camera  technique  provided  immense  flexibility  optimizing  the  intensities 
of  both  paints  simultaneously.  If  a  color  camera  was  utilized,  the  quality  of  the  PSP  and 
TSP  data  suffered  greatly  due  to  the  need  for  adequate  emission  signals  in  both  channels 
at  all  operating  conditions  while  also  not  maxing  out  the  pixel  intensity.  Additionally,  it 
was  found  that  a  typical  Bayer  RGB  color  camera  filter  did  not  adequately  separate  the 
green  and  red  emission  signals. 
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Appendix  A:  Data  Acquisition  Post  Processing  Code 


%%  Initialize  Code 
clc 

close  all 
clear  all 
tic 

for  run=2%:28; 

Chord=0 . 127 ; %  meters 
freqE=4000;%  DAQ  board  frequency 
%%  ATM 

Pa= [102400  102400  102400  102400  102400  102400  102400  102400  102400 
102400  102300  102300  102300  102300  102300  102300  102300  102300  102200 

102200  102200  102200  102200  102200  102200  102300  102300  102300  ];%  [Pa] 

Freq= [ 0  0  2.1  9.5  15.25  21  0  9.5  9.5  2.1  15.25  21  0  21  2.1  9.5 

15.25  0  9.5  2.1  15.25  21  0  2.1  9.5  15.25  21  0];  %  [Hz] 

ImageE= [ 0  100  59  59  59  59  100  59  59  59  59  59  100  59  59  59  59  100  59 
59  59  59  100  59  59  59  59  100  0000];  %  Number  of  images  per  run 
Pa=Pa . *0 . 000145037738 ;  %  Pa  to  psi 
%%  Import  Data 
Patm=Pa (run) ; 

Pressure=importdata (strcat ( ' \Run ' , int2str  (run) ,  ' Pressure . txt ' ) ) ; 
Pressure=Pressure . data+Patm; 

Encoder=importdata (strcat ( ' \Run ' , int2str (run)  ,  ' Encoder . txt ' ) ) ; 
Encoder=Encoder . data; 

le=length (Encoder ) ;  lp=length (Pressure) ; 

%%  Calculate  the  thermistor  temperatures 
Vgl=Encoder (  :  ,  4  )  ; 

Vg2=Encoder ( : , 5)  ; 

Rl=wheat stone (47000,47000,27000, 24, Vgl) ; 

T1  =  605. *R1 . A (-0 . 1266) -154 .2; 

R2=wheat stone (47000,47000,27000, 24, Vg2) ; 

T2  =  605. *R2 . A (-0 . 1266) -154 .2; 

%%  Name  the  encoder  data  and  apply  calibrations  to  the  data 
po=Encoder (1 :le,  1)  .*9. 98-. 0266;  ps=Encoder  (1 :le,  2) *5.0442-0.1; 
[mpo,sI]  =  max(po); 

ze=Encoder (1 : le,  10)  ;  ae=Encoder (1 : le,  9)  ;  Shutterl=Encoder (1 : le, 6) ; 
Trig=Encoder (1 : le, 7)  ; 

To=Encoder (1 : le, 8) ;  Shutter2=Encoder (1 : le, 3) ; 

Et= (1 : length (po) ) . / freqE; 
x— 1 : le; 

%%  low  pass  filter  pressures 

[b,a]  =  cheby2 (1, 25, . 5) ;  %  create  filter 

po  =  filtfilt (b, a, po) ;  %  filter  the  total  and  static  pressures 
ps  =  filtfilt (b, a, ps) ; 

%%  calculate  the  phase  shift  for  the  unsteady  Mach  oscillation  runs 
if  Freq(run)==0 
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else 


[b,a]  =  cheby2 (2, 35, . 06) ;  %  create  filter 
pof  =  filtfilt  (b, a,po) ;  %  use  filter  on  Po  and  Ps 
psf  =  filtfilt (b, a, ps) ; 
nnn=3 ; 

tx=-l* (ceil (nnn/Freq (run) *freqE) )  : ceil (nnn/Freq (run) *freqE)  ; 

shifter=xcov (pof (sl:sl+ceil (nnn/Freq (run) *freqE) ) , psf (sl:sl+ceil (nnn/Fr 
eq (run) *freqE) ) ) ; 

[mx, ix]  =  max ( shifter ) ; 
lag  =  tx ( ix) ; 
po=circshift (po, -lag) ; 

end 

%  add  in  the  shift  for  the  distance  between  the  Ps  and  test  section 
po=circshift (po, -7) ;  ps=circshift (ps, -7) ; 

m=real (sqrt (5* ( (po . /ps)  . A  (2/7 ) -1) ) ) ; %  calculate  Mach  number 
dm= (m-circshift (m, 1) ) *4000; %  calculate  the  rate  of  change  of  M 
%%  Calculate  the  time  stamp  for  the  pressure  brick  data 
Pt=Pressure ( 1 : lp, 1 ) ./1000; 

Po=Pressure ( 1 : lp, 2 ) ;  %  total  pressure 

Ps= ( Pressure ( 1 : Ip, 3 ) +Pressure ( 1 : lp, 4 ) ) . /2 ;  %  average  for  static 
pressure 

M=sqrt ( 5* ( ( Po . /Ps ) . A (2 /7 ) -1 ) ) ;  M=real(M);%  calulate  mach  number 
from  pressures 

%%  Calculate  the  angular  encoder  data  from  the  once  per  rev  signal 

(ze) 

%  and  500  per  rev  signal  (ae) 
dze=ze-circ shift (ze, 1) ;dze (1) =0; 
dae=ae-circ shift (ae, 1 ) ; dae ( 1 ) =0 ; 
pulse=f ind (dze) ; 
lpulse=length (pulse) ; 
angle=zeros (size  (ae) ) ; 
for  i=l : lpulse-1 ; 

angle (pulse ( i ) : pulse ( i+1 ) -1 ) =cumsum (dae (pulse ( i ) :pulse(i+l)~ 

1)  )  ; 

end 

angle=angle . *360 . /250 ; %  convert  the  signal  to  degrees 

sse=f ind (m (pulse) >0 . 4 , 1 , 'first') 

see=f ind (m (pulse) >0 . 4 , 1 ,  'last') 

mangle=0 ; 

for  i=l : see-sse-1 

[C, II ( i ) ] =max (m ( pulse ( sse+1 ) : (pulse (sse+2)- 
pulse (sse+1) ) /2+pulse (sse+1) ) )  ; 

mLi=pulse (sse+1) +11  (i) ; 
mangle (i) =angle (mLi) ; 

end 

shif ter=median (mangle) -90 
asmuth=angle-shif ter  ; 
lowasmuth= (asmuth  <  0)*360; 
highasmuth2= (asmuth  >  360)*360; 
asmuth=asmuth+lowasmuth-highasmuth; 

%%  find  camera  pulses 
cam=Shutterl>l  ; 

Cs=(cam>0)  &  (circshift (cam, 1 ) ==0 ) ; 

Ce=(cam==0)  &  (circshift (cam, 1) >0) ; 
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numb=length (find (Ce)  )  ; 
eind ( 1 : numb, 1 ) =f ind (Cs )  ; 
length ( f ind (Cs ) ) ; 
eind ( 1 : numb, 2 ) =f ind (Ce)  ; 

TimeS ( 1 : numb) =Et (eind ( 1 : numb, 1 ) ) ; 

TimeE ( 1 : numb) =Et (eind ( 1 : numb, 2 ) ) ; 
tleng=TimeE . *1000-TimeS .  *1000 ; 

Angle ( 1 : numb, 1 ) =angle (eind ( : , 1 ) ) ; 

%%  Caculate  the  velocity  from  Mach  number 
[b,a]  =  cheby2 ( 1 , 50 , . 2 ) ; 

To  =  f iltf ilt (b, a, To)  ; 

To=(To+273.1) ; %degK 
R=287 . 1; 

aa= (sqrt (1.4.*R.*To) )  ; 
vel=m. *aa; 

%%  Find  pressure  index  corisponding  to  camera  pulse  time  stamp 
for  i=l : length (TimeS) 

pind ( i , 1 ) =f ind ( Pt>TimeS ( i )  ,  1 )  ; 
pind ( i , 2 ) =f ind ( Pt>TimeE ( i )  ,  1 )  ; 

end 

Cp= zeros (size (Pressure (1 : Ip, 5 : 15) ) ) ; 
g= 1.4; 
for  i=l:ll 

Cp ( :  ,  i ) =2 . * ( Pressure (:,i+4)-Ps(:,l))  ./ (g . *M ( : , 1 )  . *Ps ( : , 1) ) ; 

end 

toe 

row=po ./0. 000145037738. /(To) ./R; 
visc=l . 458E-6 . *To . A1 . 5  .  /  (110.4+To)  ; 

Re=0 . 127 . *vel . *row. /vise; 
image=l : length (pind)  ; 
for  i=l : length (pind) 

mOut ( i ) =median (m (eind ( i , 1 ) : (eind ( i , 2 ) ) ) ) ; 
dmOut ( i ) =median (dm (eind ( i , 1 )  :  (eind ( i ,  2  )  )  )  )  ; 

ToOut (i) =median (To (eind  (i, 1)  :  (eind (i, 2) ) ) ) ; 
poOut (i) =median (po (eind ( i , 1 ) : (eind (i, 2) ) ) ) ; 
psOut (i) =median (ps (eind (i, 1) : (eind (i, 2) ) ) ) ; 

EtOut ( i ) =median (Et (eind ( i , 1 ) : (eind (i, 2) ) ) ) ; 

AngleOut (i) =median (angle (eind (i, 1) : (eind (i, 2) ) ) ) ; 

TIOut (i) =median (T1 (eind (i, 1) : (eind (i, 2) ) ) ) ; 

T20ut (i) =median (T2 (eind (i, 1) : (eind (i, 2) ) ) ) ; 

Shutter20ut (i) =median (Shutter2 (eind (i, 1) : (eind (i, 2) ) ) ) ; 
asmuthOut (i) =median (asmuth (eind (i, 1) : (eind (i, 2) ) ) ) ; 

ReOut (i) =median (Re (eind ( i , 1 ) : (eind (i, 2) ) ) ) ; 

end 

Pdata (1, :)={ 'Image',  '05',  '10',  '15',  '25',  '30',  '35',  '40',  '50',  '70',  '90',' 
-10 ' ,  ' -50 '  ,  'Mach'  ,  ' To '  ,  ' Po '  ,  ' Ps '  ,  ' Time ' ,  ' Ang' ,  ' T20  '  ,  '  T85  '  ,  '2', 

' dMach ' , ' asmuth ' } ; 

Pdata (2 : numb+1 ,:)=[ image '  Pressure (pind (:, 1 ), 5 : 1 6 )  mOut '  ToOut' 
poOut'  psOut'  EtOut'  AngleOut'  TIOut'  T20ut'  Shutter20ut'  dmOut' 
asmuthOut '  ReOut ' ] ; 

xlswrite ( (streat ( ' Run  , int2str (run) , ' PressureData ' ) ) ,  Pdata) 

PdataNames (1,  : ) ={ ' Image' ,  ' 05  ,  ' 10 '  ,  15,  25,  '30 ','35 ','40',  '50 ','70', 
90',  '-10',  '- 
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50 ' ,  ' Mach '  ,  ' To ' ,  ' Po ' ,  ' Ps ' ,  ' Time ' ,  ' Ang '  ,  ' T2  0 '  ,  ' T85 '  ,  '2',  ' dMach ' ,  ' asmuth ' 

,  '  Re  '  }  ; 

xlswrite ( (strcat ( ' Run : , int2str (run)  ,  ' PressureData ' ) )  ,  PdataNames) 
pmean= [ 1  19  20  39  40  59]; 

Pmeanl2 (1, : ) = [mean (Pressure (pind (pmean (1) , 1)  rpind (pmean (2) ,1) ,5:16) ) 

mean (m (eind (pmean ( 1 ) , 1 ) : (eind (pmean (2 ) , 2 ) ) ) ) 

mean (To (eind (pmean (1) ,1) : (eind (pmean (2) , 2) ) ) ) 

mean (po (eind (pmean (1) ,1) : (eind (pmean (2) , 2) ) ) ) 

mean (ps (eind (pmean (1) ,1) : (eind (pmean (2) , 2) ) ) ) 

mean (T1 (eind (pmean (1) ,1) : (eind (pmean (2) , 2) ) ) ) 

mean (T2 (eind (pmean (1) ,1) : (eind (pmean (2) , 2) ) )  )  ]  ; 

Pmeanl2 (2, : ) = [mean (Pressure (pind (pmean (3) , 1)  rpind (pmean (4) ,1) ,5:16) ) 

mean (m (eind (pmean ( 3 ) , 1 ) : (eind (pmean ( 4 ) , 2 ) ) ) ) 

mean (To (eind (pmean (3) ,1) : (eind (pmean (4) , 2) ) ) ) 

mean (po (eind (pmean (3) ,1) : (eind (pmean (4) , 2) ) ) ) 

mean (ps (eind (pmean (3) ,1) : (eind (pmean (4) , 2) ) ) ) 

mean (T1 (eind (pmean (3) ,1) : (eind (pmean (4) , 2) ) ) ) 

mean (T2 (eind (pmean (3) ,1) : (eind (pmean (4) , 2) ) ) ) ] ; 

Pmeanl2 ( 3 , : ) = [mean ( Pressure (pind (pmean ( 5 ) , 1 ) : pind (pmean (6) ,1) ,5:16) ) 

mean (m (eind (pmean ( 5 ) , 1 ) : (eind (pmean ( 6 ) , 2 ) ) ) ) 

mean (To (eind (pmean (5) ,1) : (eind (pmean (6) , 2) ) ) ) 

mean (po (eind (pmean (5) ,1) : (eind (pmean (6) , 2) ) ) ) 

mean (ps (eind (pmean (5) ,1) : (eind (pmean (6) , 2) ) ) ) 

mean (T1 (eind (pmean (5) ,1) : (eind (pmean (6) , 2) ) ) ) 

mean (T2 (eind (pmean (5) ,1) : (eind (pmean (6) , 2) ) ) ) ] ; 

save (['Run'  int2str(run)  ' \Pmeanl2 .mat ' ] ,  'Pmeanl2'); 
clear  all 

end 
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Appendix  B:  Image  Registration  Code 


%%  Start 
clc 

close  all 
clear  all% 
tic 

%%  for  looping 
for  dateN=  1 : 2 

%%  pick  a  date 
if  dateN==l 

date=' 11 . 21 . 11 ' ;  runS=2;  runE=28; 

else 

date=  11.23.11';  runS=l;  runE=32; 

end 

%%  for  looping 
for  run=runS : runE 

%%  Number  of  images  per  run 

ImageE= [ 0  100  59  59  59  59  100  59  59  59  59  59  100  59  59  59  59 
100  59  59  59  59  100  59  59  59  59  100  0  0  0  0;100  59  59  59  59  59  100  59 
59  59  59  100  500  500  500  500  500  500  50  50  50  50  50  50  50  50  50  50  50 
50  50  50]; 

disp([date  '  Run  '  int2str (run) ] ) 

%%  for  looping 

for  image=l : ImageE (dateN, run) ; % 

%%  initialize 
Buf fer=6 ; 
load  PSPcolormap; 
f ontsize  =  13 ; 
if  date='  11 .21 . 11 ' 

A= [ 9  9  9  9  9  9  4.5  4.5  10  10  10  10  10  0  0  0  0  0  -9  -9  - 

9  -9  -9  -10  -10  -10  -10  -10]  ; 

Freq= [ 0  0  2.1  9.5  15.25  21  0  9.5  9.5  2.1  15.25  21  0  21 
2.1  9.5  15.25  0  9.5  2.1  15.25  21  0  2.1  9.5  15.25  21  0]; 

Pa= [102400  102400  102400  102400  102400  102400  102400 
102400  102400  102400  102300  102300  102300  102300  102300  102300  102300 
102300  102200  102200  102200  102200  102200  102200  102200  102300  102300 
102300] ; %Pa 

Ta= [7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.8  7.8  7.8 

7.8  7.8  7.8  7.8  7.8  7.8  7.8  7.8  7  7  7  7  7.2  7.2  7.2];%deg  C 

Pa=Pa .*0.000145037738;  %Pa  to  psi 

else 

A=  [-4.5  -4.5  -11  -11  -11  -11  -11  11  11  11  11  11  11  11 

10  10  9  9  9  10  11  -11  -10  -9  0  0  9  10  11  -11  -10  -9] ; 

Freq= [ 0  9.5  2.1  9.5  15.25  21  0  2.1  9.5  15.25  21  0  2  9.5 
9.5229.500000000000000]; 
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Pa= [101800  101800  101800  101800  101900  101900  101900 
101900  101900  101900  101900  101900  101900  101900  101900  101900  101900 
101900  102000  102000  102000  102000  102000  102000  102000  102000  102000 
102000  102000  102000  102000  102100] ;%Pa 

Ta=  [ 5  5  5  5  6.1  6.1  6.1  6.1  6.7  6.7  6.7  8  8  8  8  8  8  9.4 
9.4  9.4  9.4  10  10  10  10  10  10  10  10  10  10  8.9];%deg  C 
Pa=Pa .*0.000145037738;  %Pa  to  psi 

end 

%%  Importing  the  pressure  data 
Ai=A (run) ; 


data=importdata ( [date, ' \Run ' , int2str (run) , ' PressureData . xls ' ] ) ; 
Pdata=data . data . Sheetl ; 

%%  Import  images 
if  image<10 

darkG  =  double (imread ( [date, ' \Run ' , int2str (run) , 
' \Dark_000 ' , int2str (image)  ,  ;  . tif ' ] ) ) ; 

onG  =  double (imread ( [date, ' \Run 1 , int2str (run) , 
' \WindOn_000 ' , int2str (image) ,  '  . tif ' ] ) )  ; 

offG  =  double (imread ( [date, ' \Run 1 , int2str (run) , 
' \WindOff_000 ' , int2str (image) ,  '  . tif '  ]  )  )  ; 

darkR  =  double (imread ( [date, ' \Run  , int2str (run) , 
' Red\Dark_000 1 , int2str (image)  ,  '  . tif ' ] ) )  ; 

onR  =  double (imread ( [date, ' \Run ! , int2str (run) , 
' Red\WindOn_000 ' , int2str (image) ,  '  . tif  '  ]  )  )  ; 
of  f  R  = 


double (imread ( [date, ' \Run ’ , int2str (run) , ' Red\WindOff  000 ' , int2str (image 

) ,  '  .tif ' ] ) ) ; 


elseif  image<100 

darkG  =  double (imread ( [date, 
' \Dark_00 ' , int2str (image) , ' . tif ' ] ) ) ; 

onG  =  double (imread ( [date, 
' \WindOn_00 ' , int2str (image)  ,  '  . tif ' ] ) ) ; 

offG  =  double (imread ( [date, 
' \WindOff_00 ' , int2str (image) , 1 . tif ' ] ) ) ; 

darkR  =  double (imread ( [date, 
' Red\Dark_00 ' , int2str (image)  ,  '  . tif ' ] ) ) ; 

onR  =  double (imread ( [date, 
' Red\WindOn_00 ' , int2str (image) , ' . tif ' ] ) ) ; 
of  fR  = 


' \Run ' , int2str (run) , 
' \Run ' , int2str (run) , 
' \Run ! , int2str (run)  , 
' \Run ' , int2str (run)  , 
' \Run ' , int2str (run) , 


double (imread ( [date,  ' \Run '  ,  int2str (run)  ,  Red\WindOf f  00 ' , int2str (image) 

, '.tif'])); 


else 


darkG  =  double (imread ( [date, 
' \Dark_0 ' , int2str (image)  ,  1  . tif ' ] ) ) ; 

onG  =  double (imread ( [date, 
' \WindOn_0 ' , int2str (image)  ,  ' . tif ' ] ) ) ; 

offG  =  double (imread ( [date, 
' \WindOff_0 ' , int2str (image)  ,  '  . tif ' ] ) ) ; 

darkR  =  double (imread ( [date, 
' Red\Dark_0 ' , int2str (image)  ,  ' . tif ' ] ) ) ; 

onR  =  double (imread ( [date, 
' Red\WindOn  0 ; , int2str (image) , ' . tif ' ] ) ) ; 


' \Run ! , int2str (run) , 
' \Run ! , int2str (run)  , 
' \Run ; , int2str (run)  , 
' \Run ' , int2str (run) , 
' \Run  , int2str (run) , 


110 


of  f  R  = 

double (imread ( [date, ' \Run  , int2str (run) , ' Red\WindOff  0 ' , int2str (image) 

'  .tif '  ] ) )  ; 

end 

%%  dark  correction  and  Resize  to  600  x  800  pixels 

onR= (onR-darkR) ; 

of fR= (of fR-darkR) ; clear  darkR; 

onG=f liplr (onG-darkG) ; 

of fG=f liplr (of fG-darkG) ; clear  darkG; 

onR  =  imresize (onR,  [600  800],  'bilinear'); 

offR  =  imresize (of fR,  [600  800],  'bilinear'); 

onG  =  imresize (onG,  [600  800],  'bilinear'); 

offG  =  imresize (of fG,  [600  800],  'bilinear'); 

[m,  n]  =  size (onG); 

%%  Std  and  Image  Rotation 
onRStd=medf ilt2 (single (stdfilt (onR)  )  )  ; 
of f RStd=medf ilt2 (single (stdfilt (offR) ) ) ; 
onGStd=medf ilt2 (single (stdfilt (onG) ) )  ; 
of fGStd=medf ilt2 (single (stdfilt (offG)  )  )  ; 

SNRr=offR. / stdfilt (offR)  ; 

SNRg=offG. / stdfilt (offG) ; 
maskOnG  =  zeros (m,  n) ; 

maskOnG(:,:)  =  single (onGStd (:,: )  <  max (max (onGStd) ) /4 ) ; 
maskOffG  =  zeros (m,  n) ; 
maskOf fG ( ; , : )  =  single (of fGStd (:,: )  < 
max (max (of fGStd) ) /4 ) ; 

of fGStd=of fGStd . *maskOffG; 
onGStd=onGStd . *maskOnG; 
var=onRStd; 

var=medfilt2 (var, [5  5] ) ; 

angle ( 1 ) =atan ( ( find (var ( 500 , : ) >max (var (500,1:250) ) * . 9, 1 , 'first ') - 
find (var (100, : ) . . . 

>max (var ( 100 , 1:250) ) *.9,1,  'first' ) ) / 4  0  0 ) ; 
onR=imrotate (onR, -l*angle (1) *180/pi ()  ,  'bilinear'  ,  ' crop ' )  ; 
var=of fRStd; 

var=medfilt2 (var, [5  5] ) ; 

angle (2 ) =atan ( ( find (var ( 500 , : ) >max (var(500, 1:250) ) * . 9, 1 , 'first ') - 
find (var (100, : ) . . . 

>max (var ( 100 , 1:250) ) *.9,1, 'first' ) ) /400) ; 
of f R=imrotate (offR, -1  * angle (l)*180/pi(),  'bilinear’,  1  crop ' ) 
var=onGStd; 

var=medfilt2 (var, [5  5] ) ; 

angle ( 3 ) =atan ( ( find (var ( 500 , : ) >max (var(500, 1:250) ) * . 9, 1 , 'first ') - 
find (var (100, : ) . . . 

>max (var (100,1:250) ) *.9, 1, 'first' ) ) /400) ; 
onG=imrotate (onG,  -  l*angle (3)*180/pi()  ,  'bilinear',  ' crop ' ) ; 
var=of fGStd; 

var=medfilt2 (var, [5  5] ) ; 

angle ( 4 ) =atan ( (find (var ( 500 , : ) >max (var (500, 1:250))*. 9,1, 'first ') - 
find (var ( 100 , : ) . . . 

>max (var ( 100 , 1:250) ) *.9, 1, 'first' ) ) /400) ; 
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rotation 
boarder) ) 
boarder) ) 
boarder) ) 
boarder) ) 


of fG=imrotate (of fG, -l*angle  (3)*180/pi(),  bilinear',  ' crop ' ) 

clear  onRStd 

clear  offRStd 

clear  onGStd 

clear  offGStd 

clear  SNRr 

clear  SNRg 

%%  Std  with  boader  cropping 

boarder=6;  %Std  while  removing  the  extra  data  from  the 

onRStd=double ( stdf ilt (onR (boarder :m-boarder , boarder : n- 

of f RStd=double ( stdf ilt (of f R (boarder :m-boarder , boarder : n- 

onGStd=double ( stdf ilt (onG (boarder :m-boarder, boarder : n- 

of fGStd=double ( stdf ilt (of fG (boarder :m-boarder, boarder : n- 

of f GStd=of fGStd (1:530,:); 
onGStd=onGStd (1:530,  : )  ; 
of f RStd=of f RStd (1:530,:); 
onRStd=onRStd (1:530,  : )  ; 

%%  Image  Registration  Translation 
if  image==l 

reflmage=  offRStd; 

end 


%  Red  off  registration 

[registration  Shift]  =  dftregistration (fft2 (ref Image) , 
fft2 (onRStd) ,  1) ; 

shift(l:2)  =  registration ( 3 : 4 ) ; 
onR  =  circshift (onR, shift); 


%  Red  off  registration 

[registration  Shift]  =  dftregistration (fft2 (ref Image) , 
fft2 (offRStd) ,  1); 

shift(l:2)  =  registration ( 3 : 4 ) ; 
offR  =  circshift (of fR, shift); 


%  Green  On  registration 

[registration  Shift]  =  dftregistration (fft2 (ref Image) , 
fft2 (onGStd) ,  1) ; 

shift(l:2)  =  registration ( 3 : 4 ) ; 
onG  =  circshift (onG, shift) ; 

%  Green  Off  registration 

[registration  Shift]  =  dftregistration (fft2 (ref Image) , 
fft2 (offGStd) ,  1); 

shift(l:2)  =  registration ( 3 : 4 ) ; 

offG  =  circshift (of fG, shift) ; 

clear  onRStd 

clear  offRStd 

clear  onGStd 

clear  offGStd 

%%  RatioG  +  R 

RatioR=onR. / of fR; 
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RatioG= (onG. /of fG) ; 

%%  Masking  the  background  of  the  image 
mask  =  zeros (m,  n) ; 

mask):, :)  =  single (onR )  >  max (max (onR) ) * . 08) ; 
mask  =  imresize (mask,  [600  800],  'bilinear'); 

%%  Find  leading  and  trailing  edges 
lead  =  find (mask (250 ,:)>. 5 , 1 ,' first ') ; 
trail  =  find (mask (250, :)>. 5, 1, ' last ') ; 

%%  Boarder  Cropping  using  Buffer 
RatioR=RatioR (Buffer+15 :m-Buf fer,  lead: trail)  ; 
RatioG=RatioG (Buffer+15 :m-Buf fer, lead: trail) ; 

%%  file  saving 

save ( [date, ' \Run '  int2str(run)  '\ Image'  int2str (image) 
' RatioR.mat ' ] ,  'RatioRmf'); 

save ( [date, ' \Run '  int2str(run)  '\ Image'  int2str (image) 
' RatioG.mat ' ]  ,  'RatioGmf'); 
end 
toe 

end 

end 
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Appendix  C:  Averaging  Ratio  of  Red  and  Green  Code 


o,  o. 
o  o 


Start 


clc 

clear  all% 
tic 

for  dateN=l 

if  dateN==l 

%%  for  looping 

date= ' 11 . 21 . 11 ' ;  runS=2;  runE=28; 

ImageE= [ 0  100  59  59  59  59  100  59  59  59  59  59  100  59  59  59  59 
100  59  59  59  59  100  59  59  59  59  100  0  0  0  0]; 
else 


%%  for  looping 

date= ' 11 . 23 . 11 ' ;  runS=19;  runE=32; 

ImageE= [100  59  59  59  59  59  100  59  59  59  59  100  500  500  500  500 
500  500  50  50  50  50  50  50  50  50  50  50  50  50  50  50]; 
end 

for  run=runS : runE 

MeanRatioRl  =  zeros (574, 514)  ; 

MeanRatioGl=zeros (574, 514) ; 

MeanRatioR2=zeros (574, 514) ; 

MeanRatioG2  =  zeros (574, 514)  ; 

MeanRatioR3=zeros (574, 514)  ; 

MeanRatioG3=zeros (574, 514)  ; 

Mask=ones (574, 514) ; 
for  image=l : ImageE (run) ; 

load ( [date, ' \Run '  int2str(run)  '\Image'  int2str (image) 

' RatioR.mat ' ] ) ; 

load ( [date, ' \Run '  int2str(run)  '\Image'  int2str (image) 

' RatioG.mat ' ] ) ; 


RatioR  =  imresize (RatioR,  [574  514],  'bilinear'); 
RatioG  =  imresize (RatioG,  [574  514],  'bilinear'); 
RatioGg=  medfilt2 (RatioG,  [50  30]); 

GI therm ( 1 , image) =RatioGg (564,111)  ; 

GI therm (2 , image) =RatioGg ( 553 , 445 )  ; 
i=0;  j=0;  k=0; 
if  image<20 


i=i+l ; 

MeanRatioRl=MeanRatioRl+RatioR; 
MeanRatioGl=MeanRatioGl+RatioG; 
elseif  image<40 

j=j+i; 


MeanRatioR2=MeanRatioR2+RatioR; 
MeanRatioG2=MeanRatioG2+RatioG; 
elseif  image<60 
k=k+l ; 
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MeanRatioR3=MeanRatioR3+RatioR; 

MeanRatioG3=MeanRatioG3+RatioG; 


else 

end 

RatioRstd=medf ilt2 (single (stdfilt (RatioR) ) )  ; 
SNRr=RatioR . /RatioRstd; 

SNRr=  medfilt2 (SNRr,  [45  30] )  ; 

MaskI=SNRr>25; 

Mask=Mask . *MaskI ; 


end 

MeanRatioRl 
MeanRatioGl 
MeanRatioR2 
MeanRatioG2 
MeanRatioR3 
MeanRatioG3 
save ( [date, 
save ( [date, 
'MeanRatioRl ' ) ; 

save ( [date, 

' MeanRatioGl ' ) ; 

save ( [date, 
'MeanRatioR2 ' ) ; 

save ( [date, 
'MeanRatioG2 ' ) ; 

save ( [date, 

' MeanRatioR3 ' ) ; 

save ( [date, 

' MeanRatioG3 ' ) ; 

save ( [date, 
toe 


:MeanRatioRl .  /  i  ; 
:MeanRatioGl . /i  ; 
:MeanRatioR2  .  /  j  ; 
:MeanRatioG2  .  /  j  ; 
:MeanRatioR3  .  /k; 
:MeanRatioG3  .  /k; 


\Run  ' 
\Run  ' 

int2str (run 
int2str (run 

\Run  ' 

int2str (run 

\Run  ' 

int2str (run 

\Run  ' 

int2str (run 

\Run  ' 

int2str (run 

\Run  ' 

int2str (run 

\Run  ' 

int2str (run 

' \Mask.mat ' ] ,  'Mask'); 

' \MeanRatioRl . mat ' ] , 

' \MeanRatioGl . mat ' ] , 

' \MeanRatioR2 . mat ' ]  , 

' \MeanRatioG2 . mat '  ]  , 

' \MeanRatioR3 . mat '  ]  , 

' \MeanRatioG3 . mat ' ] , 

' \GItherm . mat ' ] ,  'GItherm') 


end 
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Appendix  D:  In  situ  Calibration  and  Implementation  Code 


o,  o, 
o  o 


Start 


clc 

close  all 
clear  all% 
tic 

font  size=14; 
line  width=2; 

set ( 0 , ' Def aultAxesBox ' , ' on ' , ' Def aultAxesLineWidth ' , line  width, ' Default A 
xesFontSize ' , font_size) ; 

set ( 0  ,  ' Def aultLineLineWidth ' , line  width) ; 

%%  for  loop  that 
for  dateN=l 

if  dateN==l 


%%  for  looping 

date=' 11 . 21 . 11 ' ;  runS=2;  runE=28; 

Pa= [102400  102400  102400  102400  102400  102400  102400  102400 
102400  102400  102300  102300  102300  102300  102300  102300  102300  102300 
102200  102200  102200  102200  102200  102200  102200  102300  102300  102300 


]  ;%Pa 

Ta= [7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.8  7.8  7.8  7.8  7.8 
7.8  7.8  7.8  7.8  7.8  7.8  7  7  7  7  7.2  7.2  7.2];%deg  C 
Pa=Pa . *0 . 000145037738 ;  %Pa  to  psi 

A= [ 9  9  9  9  9  9  4.5  4.5  10  10  10  10  10  0  0  0  0  0  -9  -9  -9  -9  -9 
-10  -10  -10  -10  -10]; 
else 


%%  for  looping 

date= ' 11 . 23 . 11  '  ;  runE=32;  runs=[24]; 


Pa= [101800  101800  101800  101800  101900  101900  101900  101900 


101900  101900  101900  101900 
102000  102000  102000  102000 
102000  102000  102000  102100; 

Ta= [ 5  5  5  5  6.1  6.1 


101900 
102000 
;  %Pa 
6.1  6.1 
10  10  8 


101900 

102000 


101900 

102000 


101900 

102000 


6.7  6.7 

9.4  10  10  10  10  10  10  10  10  10  10  8.9];%deg 
Pa=Pa .*0.000145037738;  %Pa  to  psi 
A= [-4.5  -4.5  -11  -11  -11  -11  -11  11 
9  10  11  -11  -10  -9  0  0  9  10  11  -11  -10  -9]; 


6.7 

C 


101900 

102000 


101900 

102000 


9.4  9.4  9.4 


11  11  11  11  11  11  10  10  9  9 


end 


for  run=runs ; %runS : runE 

%%  Pressure  Data  importing  and  naming 

ImageE= [ 0  100  59  59  59  59  100  59  59  59  59  59  100  59  59  59  59 
100  59  59  59  59  100  59  59  59  59  100  0  0  0  0;100  59  59  59  59  59  100  59 
59  59  59  100  500  500  500  500  500  500  50  50  50  50  50  50  50  50  50  50  50 
50  50  50]  ; 

image=l : ImageE (dateN, run) ; 

data=importdata ( [date,  \Run ' , int2str (run) , ' PressureData . xls ' ] ) ; 
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Pdata=data . data . Sheetl ; 

M  =Pdata (image+1 , 14 ) ; 

Patm=Pa (run) ; 
po  =Pdata ( image+1 , 1 6 ) ; 
ps  =  Pdata (image+1, 17) ; 

%%  calibration  mean  ratios 

load ( [date,  ' \Run '  int2str(run)  ' \MeanRatioRl . mat ' ] ) 
load ( [date,  ' \Run '  int2str(run)  ' \MeanRatioGl . mat ' ] ) ; 
MeanRatioGl  =  imresize (MeanRatioGl ,  [574  514],  'bilinear') 

MeanRatioRl  =  imresize (MeanRatioRl ,  [574  514],  'bilinear') 

MeanRatioRl=  medfilt2 (MeanRatioRl , [50  30]); 

MeanRatioGl=  medfilt2 (MeanRatioGl , [50  30]); 

load ( [date,  ' \Run '  int2str(run)  ' \MeanRatioR2 . mat ' ] ) 
load ( [date,  ' \Run '  int2str(run)  ' \MeanRatioG2 . mat ' ] ) ; 
MeanRatioG2  =  imresize (MeanRatioG2 ,  [574  514],  'bilinear') 

MeanRatioR2  =  imresize (MeanRatioR2 ,  [574  514],  'bilinear') 

MeanRatioR2=  medf ilt2 (MeanRatioR2  ,  [50  30]); 

MeanRatioG2=  medf ilt2 (MeanRatioG2  ,  [50  30]); 

load ( [date,  ' \Run '  int2str(run)  ' \MeanRatioR3 . mat ' ] ) 
load ( [date,  ' \Run '  int2str(run)  ' \MeanRatioG3 . mat ' ] )  ; 
MeanRatioG3  =  imresize (MeanRatioG3,  [574  514],  'bilinear') 
MeanRatioR3  =  imresize (MeanRatioR3,  [574  514],  'bilinear') 
MeanRatioR3=  medf ilt2 (MeanRatioR3 ,  [50  30]); 

MeanRatioG3=  medf ilt2 (MeanRatioG3 ,  [50  30]); 

[m  n] =size (MeanRatioG3 ) ; 

%%  load  in  the  mean  pressures 

load ([date  '\Run'  int2str(run)  ' \Pmeanl2 .mat ' ] ) % ' Pmean ' 

ptM=Pmeanl2 (1:3,1:12) ; 

mM=Pmeanl2 (1:3,13)  ; 

toM=Pmeanl2  (1:3,14)  ; 

poM=Pmeanl2 (1:3,15) ; 

psM=Pmeanl2 (1:3,16) ; 

load ( [date, ' \Run '  int2str(run)  ' \Mask.mat ' ] ) 
disp([date  '  Run  '  int2str (run) ] ) 
load  PSPcolormap; 
xoc=(l:n) ,/n.*100; 

%%  x  and  I  locations 


X  (1) 

= 

find (xoc>5 , 1 , ' 

first' ) ; 

x  (2 ) 

= 

find (xoc>10 , 1 , 

' first ' ) 

x  (3) 

= 

find (xoc>15, 1, 

' first ' ) 

x  ( 4 ) 

= 

find (xoc>25, 1, 

' first '  ) 

x  (5) 

= 

find (xoc>30 , 1 , 

' first ' ) 

x  (6) 

= 

find (xoc>35 , 1 , 

' first' ) 

x  ( 7 ) 

= 

find (xoc>40 , 1 , 

' first ' ) 

x  ( 8 ) 

= 

find (xoc>50 , 1 , 

' first '  ) 

x  ( 9) 

= 

find (xoc>7  0,1, 

' first ' ) 

x  (10 

)  = 

=  find (xoc>90 , 1 

, ' first ' 

z  (  [1 

:  3 

5  8:10] ) =338; 

z  (4) 

=283; 

z  (6) 

=390; 
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z (7) =444; 
for  i=l : 10 

I (i)  =  mean (MeanRatioRl ( z ( i ) , x ( i ) ) ) ; 

IT ( i ) =MeanRatioGl (  z  ( i )  ,  x  ( i )  )  ; 

I(i+10)  =  mean (MeanRatioR2 ( z ( i ) , x ( i ) ) ) ; 

IT (i+10) =  MeanRatioG2 (z (i)  ,  x  (i) )  ; 

I(i+20)  =  mean (MeanRatioR3 ( z ( i ) , x ( i ) ) ) ; 

IT (i+20) =MeanRatioG3 (z (i) , x (i) ) ; 

Iptap (i) =ptM (1, i) ; 

Iptap ( i+10 ) =ptM (2 , i ) ; 

Iptap ( i+20 ) =ptM ( 3 , i ) ; 

end 

%%  fitting 

ft  =  fittype (  'polyl2'  ); 
opts  =  fitoptions  (  ft  ) ; 

[fit25,  gof25]  =  fit  (  [IT',  I'],  Iptap',  ft,  opts  ); 

gof25 
f  it25 

for  image=[31  44] %1 : ImageE (dateN, run) ; 

load ( [date, ' \Run '  int2str(run)  '\Image'  int2str (image) 

' RatioR.mat ']); 

load ( [date, ' \Run '  int2str(run)  '\ Image'  int2str (image) 

' RatioG.mat ' ] ) ; 

RatioR  =  imresize (RatioR,  [574  514],  'bilinear'); 

RatioG  =  imresize (RatioG,  [574  514],  'bilinear'); 
RatioGmf=  medfilt2 (RatioG,  [50  30]); 

RatioRmf=  medfilt2 (RatioR,  [ 50  10]); 

%%  full  insitu  calibration 
Pressure=f it25 (RatioGmf , RatioRmf )  ; 

Cp= ( Pressure . *Mask . /ps ( image) -1 ) ./ (0 . 5 . *1 . 4 . *M (image) .A2) 

xu=[0  xoc (8:495) ]; 

for  i=l : length (Cp ( 3 ,:) ) 

Cp (i) =mean (Cp (325:345, i) ) ; 

end 

yu=[-2  Cp 

(8:495) ]  ; 

CpData (1, : , image) =xu; 

CpData (2 , : , image) =yu; 

pdatax= [ 5  10  15  25  30  35  40  50  70  90]; 

CpPt= (pt (image, 1:10) . /ps (image) - 
1) ./ (0.5.*1.4.*M( image) . A2 ) ; % . *Mask; 

%%  plotting 
figure 
hold  on 

plot (xu, yu, ' b ' ) 

plot (pda tax, CpPt (1 : 10) , 'bs ' ) 

set (gca, ' YDir ' , ' rev ' ) 

xlabel (' Percent  x/c') 

ylabel ( ' Cp ' ) 

axis  (  [  0  100  -3.75  1.1]) 

grid  on 

hold  off 

toe 

save ( [date, ' \Run '  int2str(run)  '\Image'  int2str (image) 

' Cp .mat ' ] ,  ' Cp ' ) ; 
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end 

save ( [date,  ' \Run '  int2str(run)  ' CpData . mat ' ] ,  'CpData') 
toe 

end 

end 
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Appendix  E:  Shock  Location  Capture  Code 


%%  Start 
clc 

close  all 
clear  all% 
tic 

for  dateN=l : 2 

if  dateN==l 

%%  for  looping 

date=  11.21.11';  runS=14;  runE=28; 

Pa= [102400  102400  102400  102400  102400  102400  102400  102400 
102400  102400  102300  102300  102300  102300  102300  102300  102300  102300 

102200  102200  102200  102200  102200  102200  102200  102300  102300  102300 

]  ;%Pa 

Ta= [7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.2  7.8  7.8  7.8  7.8  7.8 
7.8  7.8  7.8  7.8  7.8  7.8  7  7  7  7  7.2  7.2  7.2]; %deg  C 
Pa=Pa . *0 . 000145037738 ;  %Pa  to  psi 

else 

%%  for  looping 

date=  11.23.11';  runS=l;  runE=32; 

Pa= [101800  101800  101800  101800  101900  101900  101900  101900 
101900  101900  101900  101900  101900  101900  101900  101900  101900  101900 

102000  102000  102000  102000  102000  102000  102000  102000  102000  102000 

102000  102000  102000  102100] ;%Pa 

Ta=  [ 5  5  5  5  6.1  6.1  6.1  6.1  6.7  6.7  6.7  8  8  8  8  8  8  9.4  9.4  9.4 
9.4  10  10  10  10  10  10  10  10  10  10  8.9];%deg  C 
Pa=Pa . *0 . 000145037738 ;  %Pa  to  psi 

end 

for  run=runS : runE 

ImageE= [ 0  100  59  59  59  59  100  59  59  59  59  59  100  59  59  59  59 
100  59  59  59  59  100  59  59  59  59  100  0  0  0  0;100  59  59  59  59  59  100  59 
59  59  59  100  500  500  500  500  500  500  50  50  50  50  50  50  50  50  50  50  50 
50  50  50]; 

disp([date  '  Run  '  int2str (run) ] ) 
for  image=l : ImageE (dateN, run) ; 

load ( [date, ' \Run '  int2str(run)  '\ Image'  int2str (image) 

' RatioG.mat ']); 

load ( [date, ' \Run '  int2str(run)  '\ Image'  int2str (image) 

' RatioR.mat ' ] ) ; 

RatioGmf=  medfilt2 (RatioG, [40  10]); 

RatioRmf=  medfilt2 (RatioR,  [40  10]); 

RatioGcp=RatioG ( 335 ,:); clear  RatioG 
RatioRcp=RatioR ( 335 , : ) ; clear  RatioR 
RatioGcpMf=RatioGmf ( 335 , : ) ; clear  RatioGmf 
RatioRcpMf=RatioRmf (335, : ) ; 

RatioRcpMfs (1, : ) =RatioRmf (100, : ) ; 
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RatioRcpMf s (2 ,  : ) =RatioRmf ( 150 ,  : )  ; 

RatioRcpMfs (3,  : ) =RatioRmf (200,  : )  ; 

RatioRcpMfs (4,  : ) =RatioRmf (250,  : )  ; 

RatioRcpMfs ( 5 ,  : ) =RatioRmf ( 350 ,  :  )  ; 

RatioRcpMfs (6, : ) =RatioRmf (400, : ) ; 

RatioRcpMfs (7,  : ) =RatioRmf (450,  : )  ; 

[n] =length (RatioGcp) ; 
xoc= (1 :n) /n*100; 

dRatioRcpMf =RatioRcpMf ( : ) -circshift (RatioRcpMf ( : )  ,  1)  ; 
dRatioGcpMf =RatioGcpMf ( : ) -circshift (RatioGcpMf ( : ) , 1) ; 
[mv  ml] =min (dRatioRcpMf (13 : floor (n/5) ) )  ; 

SL (image, 1) =xoc (ml+13) ; %Shock  Location  for  Pressure 
[mv  ml] =max (dRatioGcpMf (13 : floor (n/5) ) ) ; 

SL (image, 2) =xoc (ml+13) ; %Slock  Location  for  Temperature 
for  s— 1 : 7 

dRatioRcpMf s (s, : ) =RatioRcpMfs (s, : ) ' - 
circshift (RatioRcpMfs  (s,  : )  '  , 1) ; 

[mv  ml] =min (dRatioRcpMf s (s, 13 : floor (n/5) ) ) ; 

SL (image, s+2) =xoc (ml+13) ; 

end 

clear  dRatioRcpMfs  RatioRcpMfs  RatioRcpMf  RatioGcpMf 

end 

xlswrite ( ( [date  ' \Run'  int2str(run)  ' \ShockLocation ' ] ) ,  SL) 
toe 

end 

end 
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Appendix  F:  Tabulated  C„,  C„,  Ct  and  Cm  Data 


Table  4:  Coefficient  of  lift,  moment  and  nonnal  force  for  the  steady  runs  with  two 
standard  deviation  shown  for  each  PSP  measurement 


a 

M 

Cn 

Ca 

Ci 

cm 

9 

0.468  ±  0.004 

0.699  ±0.051 

-0.144  ±0.019 

0.713  ±0.050 

0.049  ±  0.098 

10 

0.450  ±0.005 

1.053  ±0.076 

-0.179  ±0.027 

1.068  ±0.076 

0.027  ±0.041 

11 

0.449  ±  0.003 

0.923  ±0.090 

-0.205  ±0.024 

0.945  ±  0.091 

0.080  ±  0.048 

9 

0.5  (NASA) 

0.862 

-0.104 

0.868 

0.0209 

10 

0.5  (NASA) 

0.913 

-0.111 

0.918 

0.0244 

11 

0.5  (NASA) 

0.91 

-0.0995 

0.912 

0.0179 

9 

0.549±  0.003 

0.886  ±0.018 

-0.159  ±0.006 

0.900  ±0.019 

0.029  ±0.009 

10 

0.547  ±0.002 

0.888  ±0.032 

-0.165±  0.009 

0.903  ±0.032 

0.034  ±0.013 

11 

0.546  ±0.002 

0.946  ±  0.042 

-0.156  ±0.012 

0.959±  0.042 

0.024  ±0.022 

9 

0.609  ±  0.003 

0.872  ±0.022 

-0.133  ±0.008 

0.882  ±0.022 

0.015  ±0.012 

10 

0.607  ±0.002 

0.849  ±  0.068 

-0.126  ±0.014 

0.858  ±0.067 

0.015  ±0.032 

11 

0.606  ±0.002 

1.571  ±0.575 

-0.093  ±0.042 

1.560  ±0.557 

-0.208±0.216 

9 

0.6  (NASA) 

0.852 

-0.0702 

0.852 

0.0232 

10 

0.6  (NASA) 

0.852 

-0.0517 

0.848 

-0.0026 
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